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Whether  tropical  forests  are  sources,  sinks,  or  neutral  with  respect  to  their  carbon 
balance  with  the  atmosphere  remains  unclear.  To  address  this  issue,  estimates  of  net 
ecosystem  exchange  of  carbon  and  energy  (NEE)  were  made  for  3  years  (1998-2000) 
using  the  eddy-covariance  technique  in  a  tropical  wet  forest  in  Costa  Rica.  Mean 
daytime  NEE  was  ca.  -18  |imol  CO2  m'^  s'' (uptake)  and  mean  nighttime  NEE  4.6  |amol 
CO2  m"^  s"'  (eflQux).  However,  because  -80%  of  the  nighttime  data  in  this  forest  were 
collected  during  laminar  flow  conditions  (<  0.2  m  s"'),  nighttime  NEE  was  likely 
underestimated.  Using  an  alternative  analysis,  mean  nighttime  NEE  increased  to  6.9 
nmol  CO2  m'^  s"' .  Incident  radiation  accounted  for  ~5 1  %  of  the  variation  in  the  daytime 
fluxes,  with  temperature  and  vapor  pressure  deficit  together  accoimting  for  another 
-20%.  This  forest  was  a  slight  negative  carbon  sink  in  1998  (-0.08  to  -1.42  t  C  ha''  y"'), 
a  moderate  sink  in  1999  (-1.65  to  -3.21  t  C  ha''  y'),  and  a  strong  sink  in  2000  (-6.1  to  - 
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8.1  t  C  ha''  y"').  This  trend  is  interpreted  as  relating  to  the  dissipation  of  warm-phase  El 
Nino  effects  over  the  course  of  this  study. 

The  efifects  of  net  radiation  (R,,),  vapor  pressure  deficit  (VPD),  and  surface 
conductances  on  energy  balance  and  evapotranspiration  (ET)  were  also  determined  for 
this  forest.  Sensible  (H)  and  latent  heat  (XE)  fluxes  were  estimated  as  the  sum  of  above- 
canopy  eddy-covariance  fluxes  and  changes  in  below-canopy  heat  profiles.  Albedo  was 
-12  %  of  incident  radiation  and  did  not  differ  seasonally.  Rn  was  significantly  different 
among  years,  explaining  -79%  of  the  variation  in  each  of  the  H  and  XE  fluxes.  The 
effects  of  VPD  did  not  explain  any  additional  variation  in  heat  fluxes.  X,E  was  always 
greater  than  H  (when  Rn  exceeded  40  W  m'^).  The  dimensionless  decoupling  coefficient, 
Q  was  always  >  0.5  and  peaked  at  0.7,  suggesting  that  ET  for  this  the  forest  was 
generally  decoupled  from  physiological  controls.  There  was  better  precision  in 
estimating  X,E  flux  using  the  Priestly-Taylor  model  rather  than  using  the  more 
physiologically-based  Penman-Monteith  model.  Annual  ET  was  54-66  %  of  bulk 
precipitation  and  utilized  -  88-97  %  of  available  energy  (R„). 
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CHAPTER  1 
INTRODUCTION 


Inverse  model  calculations  based  on  atmospheric  CO2  concentrations  and 
'^C02/^^C02  and  O2/N2  ratios  indicate  that  the  terrestrial  biosphere  is  currently  a  net 
carbon  sink,  partially  offsetting  the  additions  of  CO2  from  fossil  fuel  combustion  and 
deforestation  (Schimel  et  al.  2001).  Temporal  variation  in  carbon  uptake  and  emissions 
by  terrestrial  ecosystems  has  an  effect  on  interannual  variations  in  atmospheric  CO2 
concentrations  (Schimel  et  al.  2001),  although  the  magnitude  remains  uncertain 
(Houghton  et  al.  1998,  Houghton  1996,  IPCC  1995).  Recent  attempts  to  measure  forest- 
level  carbon  exchange  with  the  atmosphere  have  focused  on  temperate,  boreal,  and  arctic 
ecosystems  (Aubinet  et  al.  2001,  Black  et  al.  2000,  Law  et  al.  2000,  Valentini  et  al.  2000, 
Clark  et  al.  1999,  Frolking  et  al  1996,  Goulden  et  al  1996,  Vourlitis  and  Oechel  1997), 
with  only  few  eddy  flux  studies  from  tropical  forests  (Fan  et  al.  1990,  Grace  et  al.  1996, 
Mahli  et  al.  1998).  However,  tropical  forests  account  for  -35%  of  global  net  primary 
productivity,  >50%  of  the  carbon  in  aboveground  terrestrial  biomass  and  -20%  of  the  soil 
carbon  (Melillo  et  al.  1993,  Dixon  et  al.  1994).  Old-growth  rainforest  were  historically 
thought  to  be  carbon  neutral  (input  =  output).  This  was  challenged  by  Grace  et  al. 
(1995a)  and  Fan  et  al.  (1990)  based  on  eddy  co variance,  who  suggested  that  apparently 
long-undisturbed  Amazon  forests  sequester  carbon.  If  generally  true,  the  implications  for 
carbon  science  and  policy  making  are  enormous. 
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The  energy  balances  of  tropical  forests  are  complex  and  many  dynamic  feedback 
mechanisms  among  radiation,  cloud  formation  and  precipitation  have  been  identified 
(Wielicki  et  al.  2002,  Hartman  et  al.  2001,  Sohn  and  Smith  1992).  This  complexity 
extends  to  the  potential  role  of  the  tropical  energy  balance  in  affecting  tropical  and  global 
climates  and  general  and  anomalous  circulations  (Kelly  and  Randall  2001,  Timmerman  et 
al.  1999,  Chen  and  Van  den  Dool  1999,  Fasullo  and  Webster  1999,  Larson  et  al.  1999). 
Much  of  our  understanding  of  these  dynamics  has  relied  on  model  results,  which  have 
shown  large  spatial  and  temporal  variability  in  both  sensible  and  latent  energy  budgets 
(KeUy  and  RandaU  2001,  Raman  et  al.  1998,  Hulme  and  Viner  1998,  Shuttleworth  1988). 
In  situ  studies  have  either  scaled  leaflevel  measurements  to  whole  canopies  (Bigelow 
2001,  Avisshar  1993,  Roberts  et  al.  1993),  or  have  estimated  the  energy  balance 
components  using  eddy  covariance  over  short  periods  (e.g.,  8  d,  Shuttleworth  et  al. 
1984).  Quantifying  the  variation  of  energy  balance  parameters  and  their  biophysical 
controls  over  longer  periods  should  allow  for  better  predictions  of  runoff  and  improved 
models  of  regional  and  global  climate. 

Both  physical  and  physiological  factors  influence  forest  energy  balance,  including 
incident  radiation,  albedo,  rain,  interception,  canopy  capacitance,  and  aerodynamic  (ga) 
and  bulk  surface  (gb)  conductances.  Incident  radiation  in  the  tropics  varies  less 
seasonally  than  to  that  at  higher  latitudes,  and  values  at  the  surface  are  more  related  to 
cloudiness  than  changes  in  solar  zenith  angle.  General  circulation  models  tend  to 
underestimate  net  radiation  in  the  tropics  because  of  uncertainties  in  estimating  surface 
albedo  and  cloud  cover  (Cramer    a/.  1 999,  Ruimy  e/ a/.  1999).  Albedos  of  forests  range 
fi-om  0. 1  -0.2,  with  annual  and  seasonal  differences  affecting  the  surface  energy  balance. 
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Large  differences  in  annual  rainfall  have  been  observed  in  the  tropics  that  are  thought  to 
be  influenced  by  el  Niiio-Southem  Oscillation  (ENSO)  and  other  anomalous  circulations. 
This  in  turn  affects  the  amounts  of  water  available  for  evapotranspiration.  A  general 
observation  is  that  -50  %  of  annual  rainfall  is  re-circulated  to  the  atmosphere  through 
transpiration  and  evaporation  of  intercepted  water,  with  the  other  50%  as  runoff 
(Shuttleworth  1989).  This  imphes  that  local  climate  is  strongly  effected  by  how  energy  is 
partitioned  at  the  surface. 

Canopy  conductances  for  tropical  forests  have  been  estimated  using  both 
ecophysiological  and  micrometeorological  approaches.  Ecophysiologists  estimate  bulk 
surface  conductance  (gb)  by  scaling  leaf-level  or  sap-flow  measurements  to  the  canopy 
(Whitehead  1998,  Dolman  et  al.  1991),  while  micrometeorologists  often  model  gb  in 
relation  to  meteorological  parameters  (Wright  et  al.  1996).  Aerodjuamic  conductance  is 
generally  calculated  as  a  function  of  horizontal  windspeed,  zeroplane  displacement,  and 
roughness  length  (Derunead  and  Bradley  1985).  Evapotranspiration  from  tropical  forests 
is  generally  thought  to  be  strongly  dependent  on  aerodynamic  conductance,  because  of 
the  high  rainfall  and  the  significant  proportion  of  the  time  when  the  canopy  is  wet, 
reducing  the  importance  of  gb  in  evaportranspiration  (Shuttleworth  1989). 

Research  objectives  are:  1)  Latent  heat  is  greater  than  sensible  heat  all  the  time 
(net  radiation  >40  W  m'^),  2)  below-canopy  energy  balance  contributes  significantly  to 
the  overall  energy  balance  across  temporal  scales,  3)  with  ecosystems  that  receive  greater 
precipitation,  the  more  total  time  the  canopy  is  wet.  Because  this  wet  tropical  forest 
receives  ~  4000  mm  y"'  of  precipitation,  more  relative  time  is  needed  to  dry  the  canopy, 
so  evapotranspiration  rates  rely  more  on  aerodynamic  conductance  than  physiological 
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controls  (bulk  canopy  conductance),  and  4)  Evapotranspiration  from  wet  tropical  forests 
is  strongly  dependant  on  net  radiation,  and  accounts  for  ~  50%  of  the  precipitation. 


CHAPTER  2 
ECOSYSTEM-LEVEL  CARBON  EXCHANGE 

The  objective  of  this  chapter  was  to  define  the  patterns  of  diurnal  and  annual  net 

carbon  dioxide  exchange  and  their  climatic  controls,  for  a  lowland  tropical  wet  forest  in 

northeastern  Costa  Rica.  At  this  forest,  large  interannual  fluctuations  in  aboveground  net 

primary  productivity  since  1984  closely  paralleled  fluctuations  in  the  atmospheric  CO2 

concentrations  and  were  strongly  negatively  correlated  with  average  nighttime  minimum 

temperatures  (Clark  et  al.  in  review).  In  this  paper,  we  assess  the  effects  of 

environmental  variations  on  the  diurnal,  seasonal,  and  interannual  patterns  of  forest-level 

carbon  exchange  for  this  forest  from  1998  through  2000. 


Materials  and  Methods 

Study  Site 

This  study  was  conducted  as  part  of  a  long  term  study  of  tropical  forest  carbon 
cycling,  the  CARBONO  project,  at  the  La  Selva  Biological  Station,  Puerto  Viejo  de 
Sarapiqui,  Costa  Rica  (10°25'  51"N,  84°  00'59"W,  elevation  80-150  m.a.s.l.).  La  Selva  is 
located  in  northeastern  Costa  Rica  in  the  Caribbean  lowlands  at  the  base  of  the  central 
volcanic  chain  and  was  classified  as  tropical  wet  forest  in  the  Holdridge  life  zone  system 
by  Hartshorn  and  Peralta  (1988).  This  forest  has  an  average  400  trees  ha"'  >  10  cm 
diameter  ha''  from  -100  species  (Lieberman  et  al.  1985),  dominated  by  the  mimosoid 
legume  canopy  species,  Pentaclethra  macroloba  (34%  of  the  basal  area,  Clark  and  Clark 


5 


6 

2000).  Mean  tree  height  is  20-25  m,  with  emergents  exceeding  60  m.  Canopy  gaps 
occupy  ~  0.01-0.04  ha  ha"'  (Denslow  and  Hartshorn  1994)  making  the  overall  canopy 
very  aerodynamically  rough.  Incident  mean  (1993-1998)  daily  solar  radiation  was  14.9 
MJ  m"^  d'',  with  a  range  from  0.4  to  31.3  MJ  m'^  d''.  Mean  annual  temperatwes  were 
24.6  °C  (1982-1998,  Organization  for  Tropical  Studies,  OTS,  unpublished  weather 
records).  Mean  annual  precipitation  is  4000  mm  (from  1963  to  2000),  with  a  short  drier 
period  from  December  to  the  end  of  May,  but  with  no  month  receiving  less  than  100  mm 
(Sanford  et  al,  1994).  Soils  range  from  relatively  fertile  Inceptisols  in  riverine  areas  to 
low  pH,  low  phosphorus  Ultisols  in  upland  areas  (Sollins  et  al.  1994). 

Moisture-laden  northeast  trade  winds  originating  over  the  Caribbean  Sea 
dominate  surface  winds  (Hassenrath  1991).  During  most  (85%)  daytime  hours,  the 
annual  mean  surface  wind  direction  is  90°.  The  wetter  season  (June  through  November) 
and  drier  season  (mid  December  through  May)  are  controlled  by  the  movement  of  the 
equatorial  low-pressure  trough,  i.e.,  the  eastern  Pacific  intertropical  convergence  zone 
(ITCZ).  During  the  drier  season,  the  sub-tropical  Hadley  cell  dominates  general 
circulations,  while  the  tropical  cell  dominates  wet  season  circulations  (Sanford  et  al. 
1994).  There  are  no  data  showing  the  exact  passing  of  the  ITCZ  in  Costa  Rica.  Hence, 
seasons  are  defined  here  as  wet  season  beginning  on  May  1  (DOY  121),  and  dry  season 
on  December  20  (DOY  353).  These  dates  are  10  d  after  the  average  date  that  the  ITZC 
passes  through  Barro  Colorado  Island,  Panama  (pers.  comm.  Steve  Paton).  Other 
circulations  may  influence  wet  season  climate  including  temporales,  polar  air  masses  that 
move  down  the  North  American  continent  generating  depressions  and  prolonged  rain 
events,  chiefly  occurring  in  November  and  December  (Schlutz  et  al.  1998).  Veranillos, 
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temporary  and  often  irregular  movement  of  the  South  Pacific  anticyclone  northward, 
create  short  dry  periods,  typicaUy  lasting  7  to  10  days  in  September  or  October.  Sanford 
et  al.  (1994)  and  Holdridge  et  a/.  (1971)  provide  fiirther  site  information  for  La  Selva, 
and  Waylen  et  al.  (1996a)  and  Hassenrath  (1991)  provide  more  details  on  its  climatology. 

Because  La  Selva  is  located  at  10°  N  latitude,  there  is  little  diurnal  change  in 
sunlight  over  the  course  of  the  year,  with  only  a  40-min  difference  in  day  length  betw^een 
solstices.  For  this  study,  sunrise  and  sunset  were  defined  as  0600  and  1 800  h,  delineating 
daytime  and  nighttime  periods. 

A  42  m  tower  (Upright  Inc.  Selma,  CA)  was  used  to  access  the  canopy 
environment  and  to  support  meteorological  instrumentation.  The  site  was  a  relatively  flat 
ridgetop  in  an  area  of  generally  rolling  topography,  with  -20-30  m  relief  between  stream 
bottoms  and  ridgetops  (OTS  unpublished  digital  elevation  model).  After  accounting  for 
stability  effects,  a  source  area  model  (Schuepp  et  al.  1 990)  was  used  to  estimate  that 
under  stable  conditions,  95%  of  the  cumulative  flux  was  derived  fi-om  within  1.2  km  of 
the  tower  (at  a  mean  horizontal  windspeed  of  3  m  s'').  The  tower  was  sited  to  minimize 
edge  effects,  below-canopy  advection  either  to  or  fi"om  the  site,  and  any  major  directional 
differences  due  to  forest  composition  and  structure. 

Meteorological  Data 

Microclimate  data  were  collected  continuously  at  the  tower  top.  Measured 
variables  included  incident  radiation  (LI- 190,  Ll-Cor  Inc.,  Lincoln,  NE),  photosynthetic 
photon  flux  density  (PPFD,  LI-200X,  Ll-Cor  Inc.),  aspirated  air  temperature  (Ta,  100  Q 
platinum  RTD,  Omega  Engineering,  Stamford,  CT),  and  bulk  rainfall  (TE525  metric. 
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Texas  Electronics,  Dallas,  TX).  Atmospheric  pressure  (PB105,  Vaisala,  Helsinki, 
Finland)  was  monitored  at  ground  level. 

All  of  the  above  data  were  collected  at  an  interval  of  5  sec  and  compiled  as  30- 
min  averages  with  dataloggers  (CRIOX  and  2 IX,  Campbell  Scientific  Inc.,  Logan,  UT). 
Instruments  were  cleaned,  leveled  as  necessary,  and  recalibrated  according  to 
manufacturers'  instructions.  At  times  when  the  PPFD  sensor  was  inoperational,  PPFD 

•     •  •    •  2 

was  estimated  by  linear  regression  equation  relating  PPFD  to  incident  radiation  (R  > 
0.99).  Likewise,  when  either  power  outages  occurred  or  aspirated  air  temperature  were 
not  logged,  air  temperature  was  estimated  fi-om  a  regression  against  an  air  temperature 
sensor  CS500  (CampbeU  Scientific)  also  mounted  at  42  m  in  a  radiation  shield  (R^  > 


Long  term  meteorological  data  fi-om  La  Selva  were  used  to  examine  decadal  scale 
trends  in  microclimate  (OTS,  http://www.ots.duke.edu).  PPFD  and  air  temperature  have 
been  measured  since  1982  and  bulk  precipitation  since  1961. 

To  assess  zero-plane  displacement  (d),  or  the  mean  level  of  momentum  absorption,  four 
3-cup  anemometers  (Model  03103-L,  R.M.  Young,  Traverse  City,  MI)  were  mounted 
vertically  along  the  tower  at  35.5,  3 1 .6,  28.2,  and  25  m  above  the  ground,  d  was  then 
estimated  by  determining  the  intercept  (yo)  of  Eq.  2-1: 


where  Z  is  measurement  height  above  the  ground  (m),  u  is  the  30-min  time  average  of  the 


0.98). 


Eq.  2-1 


instantaneous  measurement  of  horizontal  wind  velocity  (cm  s" )  at  each  height. 


Roughness  length  (zo)  was  estimated  using  d: 


(      M  *  V'  Eq.  2-2 


exp^= — 


where  Zm  is  the  measurement  height  (42  m),  u*  is  friction  velocity  (m  s" ),  Um  is  the  30- 
min  time  average  of  the  instantaneous  measurement  of  horizontal  wind  velocity  at 
measurement  height  (m  s''),  and  k  is  the  von  Karmen  constant  (0.41,  dimensionless). 
The  ratio  of  convective  to  mechanical  production  of  turbulent  kinetic  energy  (Monin- 
Obukov  length,  L)  was  used  to  determine  atmospheric  stability  as  in  Eq.  2-3: 

^  ^  -pCpT.u'  Eq.  2-3 

gkH 

where  p  is  the  density  of  air  (kg  m"^),  Cp  is  the  specific  heat  capacity  of  air  (J  kg"'  K"'),  Ta 
is  in  Kelvin  (K),  g  is  acceleration  due  to  gravity  (m  s'^),  and  H  is  the  sensible  heat  flux 
density  (J  m"^  s'')  (Pahlow  and  Parlange  2000,  Monteith  and  Unsworth  1990,  Rosenberg 
etal.  1983). 

Net  Ecosystem  Exchange  (NEE)  Measurements 

A  closed-path  eddy  co  variance  system  was  used  to  estimate  the  portion  of  NEE  of 
CO2  contributed  by  turbulent  exchange.  Because  the  below-canopy  environment  was  not 
always  subject  to  turbulent  transfer  (i.e.,  well-mixed  conditions),  a  profile  system  was 
used  to  estimate  the  rate  of  change  of  [CO2]  below  the  canopy.  NEE  was  then  estimated 
as 

^27 


NEE  =  w'C02'+         —  ^^27  Eq.  2-4 

io.5 

where  w'  and  CO2'  are  the  deviations  of  instantaneous  values  from  a  running  mean  of 
vertical  windspeed  (m  s"')  and  molar  fraction  of  CO2  (^mol  CO2  mol''),  respectively,  and 


d[C02]. 
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Zx  is  measurement  height  (m).  The  lirst  term  is  the  30-min  time-averaged  eddy 
covariance  flux.  The  second  term  is  the  storage  flux  below  42  m.  The  convention  used  is 
that  negative  values  of  NEE  correspond  to  uptake  of  CO2  by  the  forest  from  the 
atmosphere. 

A  3-D  sonic  anemometer  (K-probe,  Applied  Technologies  Inc.,  Boulder  CO)  was 
used  to  measure  wind  velocities  in  each  polar  coordinate  (w,  v,  u)  and  sonic  temperature 
(0).  The  gas  sampling  inlet  was  mounted  on  the  sonic  anemometer  and  co-located  with 
the  top  transducer  in  the  w-axis.  Infrared  gas  analyzers  (IRGA,  model  LI-6262,  Ll-Cor 
Inc.)  were  used  to  measure  concentrations  of  CO2  and  H2O  vapor,  controlled  for  pressure 
and  temperature  at  ground  level  inside  a  climate-controlled  structure.  Flow  rates  were 
maintained  by  pumps  (KNF  Neuberger,  Trenton,  N J)  and  mass  flow  controllers  (Model 
series  200,  0-10  1pm,  Teledyne  Hastings  Inc.,  Los  Angeles  CA).  Sampled  air  flowed 
through  -60  m  of  tubing  (4.8  mm  ID  Teflon  tubing)  at  a  rate  of  8  1pm  resulting  in  a 
lagtime  of -14.2  s. 

The  NOAA  flux  software  program  (McMillen  1988),  with  a  400  second  digital 
recursive  running  mean  and  a  fixed  lagtime,  was  used  to  collect  raw  eddy  covariance  data 
files.  IRGA  voltage  outputs  were  digitized  by  a  12-bit  analog-to-digital  board. 
Covariances,  wind  and  scalar  statistics,  and  coordinate  rotations  were  calculated  in  real 
time  at  1 0  hz.  Protocols  for  accuracy,  precision,  quality  control  and  assurance  were  used 
as  defined  by  the  AmeriPlux  Science  Plan 
(http://cdiac.esd.oml.gov/programs/ameriflux/scifhtm). 

The  response  of  all  instruments  must  be  as  fast  as,  or  faster  than,  the  turbulence 
that  is  carrying  the  bulk  of  mass  and  energy.  This  process  occurs  between  frequencies  of 
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1-10  hz,  within  the  inertia!  sub-range  (Kaimal  and  Finnigan  1994).  While  the  sonic 
anemometer  operates  at  these  speeds,  the  response  time  for  the  IRGAs  is  slower,  at~8-9 
hz.  To  account  for  this  frequency  loss,  fast  Fourier  transfer  (FFT)  analyses  were  applied, 
where  the  portion  of  attenuation  shown  by  the  co-spectral  density  in  the  inertial  sub-range 
was  compared  to  the  total  spectra.  This  ratio  for  w'0'  was  considered  ideal,  since 
temperature  (0)  was  measured  by  the  sonic  anemometer  at  1 0  hz,  and  then  compared  to 
the  spectra  of  w'C02'  to  estimate  a  spectral  correction  factor,  SCf: 


where  Swx'  is  spectral  density  of  w'  and  CO2'  or  9',  n  is  the  natural  frequency,  and  w'x'  is 
the  mean  covariance  of  w'  and  CO2'  or  9'  (Meyers  and  Baldocchi  1991,  Baldocchi  and 
Meyers  1 989).  Although  we  assumed  that  the  dissipation  of  turbulent  kinetic  energy 
occurred  in  the  inertial  sub-range  for  all  the  scalars,  we  did  not  expect  that  the  Kaimal 
spectral  relationship  (a  slope  of  -  2/3  within  the  inertial  sub-range)  would  hold  for  every 
30-min  period  because  of  roughness,  differing  stabilities,  and  possible  density  driven 
flows  over  time. 

Profile  measurements  were  used  to  calculate  below-canopy  CO2  storage  dynamics 
(Eq.  2-4).  CO2  was  coUected  from  6  inlets  at  0.5,  7.3,  1 1.95,  16.55,  21.2,  and  27.6  m.  A 
datalogger  (model  2 IX,  Campbell  Scientific  Inc.)  was  used  to  operate  solenoids  that 
switched  the  flow  (~3  1pm)  from  each  inlet  through  the  IRGA  (Li-Cor  6262)  for  5  min 
during  each  30-min  period  and  to  record  the  raw  data.  Platinum  resistance  thermometers 
(100  Q  PRT,  m68.  Omega  Engineering,  Stamford,  CT)  housed  in  radiation  shields  were 


Eq.  2-5 
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co-located  with  each  inlet.  When  sampling  occurred,  the  airflow  acted  to  aspirate  the 
PRTs.  Temperature  and  humidity  profiles  were  used  to  account  for  changes  in  mass  flow 
due  to  changes  in  density  (Webb  et  al.  1980).  Below-canopy  storage  was  calculated  fi-om 
line-averaged  profile  measurements  using  Eq.  2-4;  it  was  assumed  that  this  profile  was 
similar  across  the  flux  source  area. 

Both  eddy  co variance  and  profile  measurements  began  in  April  1 998  and 
continued  through  the  end  of  December  2000.  Gaps  in  the  dataset  occurred  for  periods  of 
2-14  days  when  either  power  failure  or  instrumentation  mallunctions  occurred.  IRGAs 
were  calibrated  every  2-3  days.  Improved  precision  in  calibration  was  achieved  starting 
in  February  1999  by  plumbing  nitrogen  through  the  IRGA  reference  cell  as  a  zero 
reference.  A  model  was  used  to  relate  daytime  NEE  to  PPFD  (Ruimy  et  al.  1995): 

(^>«^max  Eq.  2-6 


NEE^^^  =Re  + 

where  NEEday  was  calculated  using  Eq.  2-4,  Pmax  is  maximum  ecosystem  CO2  uptake  rate 
(fxmol  CO2  m'^  s''),  Re  is  ecosystem  respiration  (jamol  CO2  m"^  s"'),  0  is  PPFD  (jimol  m"^ 
s"'),  and  a  is  apparent  quantum  efficiency  (5CO2/50).  To  describe  the  effects  of 
temperature  on  nighttime  NEE,  a  second  model  was  used: 

where  Ro  is  the  base  ecosystem  respiration  rate  (|amol  CO2  m'^  s"')  when  air  temperature 
is  0  °C,  T  is  temperature  (°C),  and  b  is  an  empirical  coefiBcient. 

A  general  linear  model  (SAS  v.  8.01,  Gary,  NC)  was  used  to  test  first  and  second 
order  effects  of  other  variables  on  NEE,  including  PPFD,  temperature,  VPD,  season  and 
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year.  Sigmaplot  v.  5.0  (SPSS  Inc.,  Richmond,  CA)  was  used  to  describe  these 

relationships. 

Data  Screening 

Eddy  covariance  data  were  screened  for  validity  and  removed  when  either  i)  the 
standard  deviation  of  w',  CO2',  or  9'  was  >  1 .74  times  the  mean,  ii)  rain  occurred,  iii)  30- 
min  data  collection  periods  were  incomplete,  or,  iv)  signals  from  either  the  sonic 
anemometer  or  the  IRGA  were  out-of-range.  Profile  data  were  removed  when  either  i) 
data  were  beyond  3  SD  from  the  mean,  ii)  30-min  data  collection  periods  were 
incomplete,  or  iii)  signals  from  the  IRGA  were  out-of-range. 

Leaf  Area  Estimation 

Although  the  La  Selva  forest  is  lai  gely  evergreen,  seasonal  differences  in  leaf  area 
index  (LAI)  were  expected  because  8%  of  the  tree  species  are  deciduous  in  the  dry 
season,  and  28  %  of  tree  species  produce  annual  leaf  flushes,  many  at  the  onset  of  the  wet 
season  (Frankie  et  al.  1974).  Furthermore,  many  tropical  rain  forest  tree  species  are 
facultatively  deciduous,  losing  up  to  half  of  their  leaves  during  prolonged  dry  periods 
(Richards  1996).  Photographic  estimates  of  eLAI  (estimated  LAI)  were  derived  using  the 
WINPHOT  program  (Ter  Steege  1996)  each  year  during  the  wet  and  dry  seasons  across 
18  randomly  stratified  0.5  ha  plots  (description  of  statistical  design  for  plot  layout  in 
Clark  and  Clark  2000).  Within  each  plot,  6  photographs  were  made  at  each  sampling 
date  under  diSuse  light  conditions  at  the  same  randomly  chosen  points.  Because  these 
estimates  were  derived  optically  with  no  means  of  direct  calibration,  they  should  be 
viewed  relatively. 
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Characterizing  Canopy-Level  Turbulence 

There  were  -2/3  slopes  for  the  normalized  spectra  of  wind  velocities  in  the 
inertial  sub-range  during  periods  of  both  stable  and  unstable  atmospheric  conditions, 
confirming  a  transfer  of  energy  to  the  canopy  with  shear  forces  dominating  (Figure  2-1) 
The  spectral  density  decreased  during  stable  conditions  (Figure  2- IB),  as  did  the  eddy 
covariance  flux  estimates,  but  the  general  relationships  held.  Buoyancy  forces  produced 
measurable  vertical  wind  movement  at  night,  as  indicated  by  the  positive  1:1  slopes  at 
wavenumbers  <  0. 1 .  The  observed  shift  in  the  spectral  peak  between  stable  and  imstable 
conditions  was  similar  to  that  reported  in  other  studies  (Kaimal  and  Finnigan  1994, 
Anderson  et  al.  1 986).  The  spectral  correction  factor  ranged  from  1 . 1 8  to  1 .08  and 
varied  with  u*  (Figure  2-2). 

Zero-plane  displacement  (d)  and  zo  for  momentum  also  varied  with  stability 
(Table  2-1).  During  unstable  periods  (L<  -50  m),  zo  increased  to  -2.4  m,  a  long 
roughness  length  even  for  a  forest  (Hansen  1993),  with  a  mean  level  of  drag  (d)  of -22 
m.  Aerodynamic  roughness  lengths  sharply  decreased  from  slightly  unstable  conditions 
(-50  <  L  <  -10)  to  neutral  conditions  (-10  <  L  <  10),  indicating  the  quick  formation  of 
stratified  laminar  flow  and  the  decoupling  of  the  below-canopy  environment. 

The  relationship  between  NEEday  and  u*  was  linear  (Figure  2-3A).  There  was 
decrease  in  the  NEEnight  with  u*,  although  no  threshold  was  observed.  When  NEEnight 
data  were  averaged  across  u*  intervals  of  0.025  m  s"',  the  relationship  was  strongly  linear 
below  a  u*  of  0.45  m  s''  (Figure  2-3B).  No  relationship,  however,  was  foimd  between 
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the  residuals  from  the  energy  budget  and  u*  under  any  stability  conditions,  so  that  no  u* 
threshold  could  be  determined  and  used  to  filter  data,  as  has  been  done  in  most  other 
studies  (e.g.,  Clark  et  al.  1999).  We  assumed  that  the  most  accurate  estimate  of  nighttime 
turbulent  exchange  occurred  at  larger  values  of  u*(i.e.,  >  0.4  m  s"')  and  that  conversely, 
storage  occurred  at  very  low  u*  (i.e.,  <  0.05  m-  s''),  in  spite  of  the  lack  of  obvious 
thresholds.  The  consequences  of  these  assumptions  are  discussed  below. 

Diurnal  Patterns  in  NEE 

The  diurnal  pattern  of  the  CO2  storage  flux  was  very  consistent  throughout  the 
year  (Figure  2-4 A).  The  greatest  fluxes  were  observed  in  early  morning  hours,  when 
below-canopy  CO2  that  was  respired  during  the  night  and  stored  in  the  air  column  below 
the  sonic  anemometer  was  vented  or  re-fixed  through  photosynthesis.  The  magnitude  of 
these  morning  ejections  may  be  underestimated,  because  venting  may  skew  the 
distribution  of  wind  statistics  in  the  30-min  dataset,  and  so  valid  data  could  have  been 
inadvertently  removed  during  the  screening  process. 

The  maximum  average  storage  flux  was  -5.6  j^mol  CO2  m'^  s'',  which  occurred  at 
-0800  when  the  convective  boundary  layer  was  developing,  as  shown  by  the  increasing 
vertical  windspeed  in  Figure  2-4B.  Storage  fluxes  decreased  imtil  -1400,  after  which 
only  net  effluxes  from  below  the  canopy  were  observed.  This  also  coincided  with  peaks 
in  below-canopy  temperature  (Figure  2-4B)  and  vertical  windspeed.  Storage  eflflux 
increased  to  a  peak  just  after  sunset.  The  maximum  nighttime  (before  0600)  storage 
efflux  was  2.97  fimol  CO2  m"^  s'',  with  an  average  of  1.6  ±0.13  i^mol  CO2  m'^  s"'. 
Storage  generally  decreased  throughout  the  night,  along  with  temperature  and  vertical 
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windspeed.  Vertical  windspeed  and  below-canopy  temperatures  diverged  between  1630 
and  0130,  which  indicates  horizontal  advection  of  below-canopy  CO2  ofifthe  site  may 
have  occurred  and  that  a  portion  of  the  flux  may  have  been  missed  (Mahrt  et  al  2000). 
The  maximum  NEEday  based  on  the  3-yr  mean  for  each  half-hour  was  -17.3  ±0.3  |imol 
CO2  m'^  s'',  occurred  at  ~  1 130,  and  closely  followed  the  inverse  pattern  of  PPFD  (Figure 
2-5).  Nighttime  eddy  covariance  flux  was  positive  and  fairly  constant  throughout  the 
night.  At  dawn,  this  flux  sharply  decreased  (uptake  into  the  forest)  for  -30  minutes  with 
decreases  in  Monin-Obukov  length  (Eq.  2-3). 

Environmental  Controls  on  NEE 

NEEday-  NEEday  was  negatively  correlated  with  PPFD  (Figure  2-6A)  and  had  an 
estimated  mean  maximum  of  -1 8  ±  9  ^imol  CO2  m"^  s  '.  A  linear  model  that  included 
second-order  effects  of  year,  PPFD,  and  VPD  explained  ~72%  of  the  total  variation  in 
NEEday.  PPFD  alone  accounted  for  -51%  of  the  variation  (Figure  2-6),  with  no 
significant  effect  of  season.  The  light  response  function  of  Eq.  2-6  had  a     of  0.51 . 
Residuals  from  this  function  were  weakly  related  to  temperature  and  VPD  (Figvire  2-6B 
and  2-6C). 

Since  VPD  includes  temperature  as  a  component,  VPD  and  temperature  are 
correlated,  and  since  VPD  directly  affects  stomatal  conductance  (Law  et  al.  2000, 
Landsberg  and  Gower  1997),  we  examined  the  influence  of  VPD  on  NEEday  further  by 
separating  NEEday  into  three  VPD  classes  (0-0.5,  0.51-1.00,  and  >  1.0  kPa)  and  refitting 
Eq.  2-6  for  each  year.  The  results  (Figure  2-7,  Table  2-2)  show  that  the  response  function 
is  more  linear  within  a  VPD  class.  Pmax  and  Re  in  1998  were  significantly  lower  than 
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those  found  in  1999,  while  a  in  2000  was  greater  than  in  the  other  two  years  (Table  2-2). 

2  -I 

Light  compensation  points  ranged  from  1 10  to  207  |xmol  CO2  m'  s' . 
NEEnighi.  The  mean  NEE„igh,  was  4.82  ±0.6  ^imol  CO2  m"^  s"'  (mean  +1  SD).  This  may 
be  an  underestimate  because  >80%  of  the  nighttime  turbulent  exchange  measurements 
were  made  with  u*  <  0.4  m  s"'.  Using  only  data  with  u*  >  0.4  m  s"\  the  estimate 
becomes  considerably  larger,  6.98  ^mol  CO2  m'^  s'\  with  4.83  ±0.21  and  2.15  ±0.1 1 
|j,mol  CO2  m'^  s"'  (mean  ±1  SE)  the  contributions  from  eddy  covariance  and  storage 
fluxes,  respectively.  In  addition,  the  averaged  daytime  Re  values  (from  Eq.  2-6,  Table  2- 
2)  ranged  from  5.07-6.42  |amol  CO2  m"^  s"',  supporting  use  of  the  higher  NEEnight  values 
derived  here. 

NEEnight  was  weakly  related  to  temperature  over  the  entire  sampling  period 
(Figure  2-8),  with  a  Qio  of  1 .79  (p  <  0.053). 

Estimating  Annual  NEE 

1  estimated  annual  NEE  in  two  ways.  In  the  iirst  case,  half-hourly  meteorological 
data  were  used  to  drive  Eq.  2-6  for  each  year  and  VPD  class  to  derive  an  annual  NEEday. 
Then  a  fibced  NEEnight  of  6.98  fxmol  CO2  m"^  s"'  was  subtracted.  I  refer  to  these  results  as 
NEEmodeied-  In  the  second  case,  measured  NEEday  and  NEEnight  (with  u*  >  0.4  m  s'^)  were 
used,  and  only  the  gaps  were  filled  using  the  procedure  above;  results  are  referred  to  as 
NEEgap.  There  were  marked  differences  in  the  cumulative  NEE  among  years  and  method 
of  estimation  (Figure  2-9).  Annual  NEE  was  estimated  to  be  higher  in  1999  and  2000 
than  in  1998,  with  2000  the  strongest  sink  year. 
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Discussion 

Characteristics  of  the  La  Selva  Canopy 

I  assume  that  the  eddy  covariance  measurements  from  La  Selva  were  robust  and 
represent  the  nature  of  CO2  exchange  between  the  canopy  and  the  atmosphere  on  the 
basis  of  two  requirements  generally  prescribed  for  this  method:  (1)  a  consistent  energy 
cascade  in  the  inertial  sub-range  in  the  power  spectra  under  both  stable  and  unstable 
conditions  (Kaimal  and  Finnigan  1994),  and  (2)  the  measurement  height  was  at  least  1.8 
times  d  (it  was  in  fact  at  least  8  times  zo  above  d  in  unstable  conditions  and  6  zo  above  d 
under  all  other  stability  classes,  Schmid  et  al.  2000,  Monteith  and  Unsworth  1990).  I  did 
not  expect  to  see  a  well-developed  energy  cascade  during  nighttime  conditions.  It  is 
likely  that  this  was  due  to  roughness  induced  turbulence  in  the  nighttime  flows. 

It  is  unclear  why  no  u*  threshold  in  the  eddy  covariance  data  was  observed. 
Turbulence  structures  above  a  fixed  plane  (e.g.,  above  a  vegetated  canopy)  from  other 
studies  have  generally  been  described  using  Monin-Obukov  theory  (Leclerc  et  al. 
submitted)  and  measured  over  uniform  canopies  with  short  aerodynamic  roughness.  As 
roughness  lengths  become  greater  (>  1  m),  the  effect  on  turbulence  in  the  well-mixed 
layer  and  applicability  of  Monin-Obukov  theory  become  questionable  (Nakamura  and 
Mahrt  2001,  Ayotte  et  al.  1999,  Raupach  and  Finnigan  1997).  Under  these  conditions,  u* 
becomes  homogenized  over  a  broader  range  of  turbulence  lengths  and  potentially  has  less 
explanatory  power  (Nakamura  and  Mahrt  2001).  At  La  Selva,  there  were  insuflHcient 
data  to  assess  if  the  linear  relationship  between  u*  and  NEEnight  existed  at  values  >0.45  m 
s"',  or  if  the  relation  ultimately  reached  some  asymptote.  Data  that  would  have 
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contributed  toward  developing  a  u*  filter  may  also  have  been  screened  out  by  other 
criteria. 

Environmental  Controls  on  NEE 

The  average-diurnal  pattern  in  NEE  (Figure  2-5)  is  strongly  symmetrical  around 
1 130.  On  average,  the  storage  term  contributed  33  %  to  NEEmght.  The  power  spectra  for 
nighttime  eddy  covariance  indicated  a  transfer  of  mass  and  energy,  but  the  flux  was 
relatively  small  and  the  source  distance  quite  long  (1.5-2.0  km).  The  nighttime 
envirorunent  below  the  canopy  is  quite  different  from  that  at  the  tower  top  and  subject  to 
diabatic  flows  controlled  by  changes  in  air  density  and  slope  (Mahrt  1 992),  making  the 
source  area  for  the  storage  flux  more  localized.  The  greatest  variation  in  nighttime  CO2 
profiles  occurred  during  periods  with  the  most  rapid  changes  in  Ta.  The  diurnal  patterns 
in  CO2  profiles  and  respective  storage  fluxes,  however,  were  relatively  constant 
throughout  our  measurement  period.  After  the  morning  re-assimilation  of  stored  CO2 
dissipated,  NEEday  was  clearly  dominated  by  above-canopy  fluxes. 

A  coarse  estimate  of  the  amount  of  recycled  CO2  can  be  obtained  as  the  difference 
between  integrated  NEEnight  and  the  morning  storage  flux.  This  represents  a  value  of 
32%  for  the  fraction  of  integrated  NEEnight  that  was  recycled  below  42  m,  similar  to 
isotopically-derived  estimates  of  re-synthesized  CO2  from  the  Ducke  forest  near  Manaus, 
Brazil  (39%,  Sternberg  et  al.  1997)  and  within  the  range  from  another  neotropical  forest, 
Barro  Colorado  Island,  Panama  (31-38%),  Sternberg  et  al.  1989). 

The  lack  of  a  relationship  between  30-min  averages  of  NEEnight  and  temperature 
was  likely  due  to  the  combination  of  different  factors.  First,  there  are  numerous  sources 
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of  C02  for  ecosystem  respiration,  each  with  their  own  controlling  factors  (Davidson  et  al. 
1998).  Soil  respiration  can  be  influenced  by  soil  type,  water  status  and  temperature 
(Schwendenmann  et  al  submitted),  and  foliar  respiration  may  be  influenced  by  nitrogen 
and/or  photosynthate  availability  in  the  canopy  at  the  time  of  leaf  expansion.  Second, 
there  was  only  a  small  range  annual  nighttime  temperature  (<  9  °C).  Third,  the  above- 
canopy  source  area  is  integrated  over  a  larger  area  than  that  from  the  below-canopy 
environment  (Raupach  et  al.  1992),  subjecting  the  storage  flux  to  localized  biotic  effects 
or  the  possibility  of  below-canopy  advection.  Foiuth,  as  mentioned  above,  80%  of  flux 
estimates  were  made  under  conditions  with  low  u*  (<  0.4  m  s''),  questionable  conditions 
in  this  case  for  eddy  covariance.  Finally,  carbon  may  have  been  exported  from  the 
system  as  CO2,  due  to  large  morning  ventilation  events  could  not  be  quantified  with 
confidence  (Wilson  et  al.  1998,  Mahrt  1992). 

The  relationship  between  NEEnight  and  temperature  (Figure  2-8)  based  on  longer 
averaging  intervals  is  subject  to  the  same  potential  sources  of  error  as  mentioned  above 
and  may  also  lead  to  underestimates  NEEnight-  Schwendenmann  et  al.  (submitted) 
reported  soil  respiration  rates  at  La  Selva  from  upland  soils  ranged  from  3.3  to  4.3  |amol 
CO2  m'  s"  .  Given  that  soil  respiration  accounts  for  50%  of  Re  for  a  wide  range  of 
forests  (Amthor  1994,  Ryan  1991),  this  further  suggests  that  our  higher  value  (6.98  nmol 
CO2  m'  s' )  for  NEEnight  is  more  likely  correct.  Moreover,  when  I  recalculated  annual 
NEEnight  using  the  nighttime  respiration  flinction  in  Figure  2-8,  it  also  resulted  in  a  very 
low  value  (4.3  ^mol  CO2  m"^  s''),  similar  to  that  found  using  the  NEEnight  data  across  all 
u*  values,  and  again  a  value  completely  out  of  step  with  measured  soil  respiration  and 
our  other  estimation  of  NEEnight- 
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The  potential  exists  that  below-canopy  nighttime  advective  flow  contributes  error 
in  NEEnight-  However,  for  this  to  occur  there  must  exist  specific  conditions,  such  as,  a 
strong  upslope  temperature  gradient,  a  breakdown  of  below-canopy  resistance  to  flow, 
and/or  a  net  vertical  movement  of  wind  into  the  forest.  At  La  Selva  however,  there  is 
only  a  small  change  in  temperature  over  a  large  upslope  area,  there  is  considerable 
resistance  presented  by  large  below-canopy  leaf  area  and  tree  stems,  and  nighttime 
above-canopy  vertical  and  below-canopy  horizontal  windspeeds  (data  not  shown)  were 
not  much  different  from  the  expected  accuracy  of  the  sonic  anemometer  (i.e.,  0.05  m  s"'). 
This  is  an  unresolved  issue  and  these  factors  are  concerns  for  many  tower  flux  sites 
(Massman  and  Lee  2001). 

In  the  face  of  increasing  global  temperatures  (National  Academy  of  Science 
2000),  there  is  increased  focus  on  the  role  of  temperature  in  controlling  the  carbon 
dynamics  in  the  tropics.  Kindermann  et  al.  (1996)  modeled  the  effects  of  increased 
temperature  on  carbon  stores  and  vsdth  even  small  increases  in  temperature  (~0.5  °C), 
large  effluxes  of  carbon  to  the  atmosphere  are  expected.  It  is  hypothesized  that  most  of 
this  carbon  will  be  from  the  tropics  (Trumbore  et  al.  1996).  At  La  Selva,  large  year-to- 
year  fluctuations  over  the  past  16  years  in  aboveground  biomass  increments  have  been 
negatively  correlated  with  both  the  mean  nighttime  temperature  and  variations  in  annual 
fluctuations  in  atmospheric  CO2  concentrations  (Clark  et  al.  in  review).  In  this  study,  I 
found  only  a  small  temperature  influence  on  NEEday  and  NEEnight  and  only  when  all  three 
years  of  data  were  pooled  and  hourly  averages  were  used.  I  cannot  dismiss  the  NEEday 
response  to  temperatm-e  as  entirely  due  to  VPD,  and  may  in  fact  be  partially  due  to 
photorespiration.  Interestingly,  Grace  et  al.  (1996)  and  Malhi  et  al.  (1998)  did  not  report 
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a  nighttime  temperature  response  for  Amazonian  forests  using  eddy  covariance  data.  It 
may  well  be  that  effects  of  temperature  on  ecosystem  respiration  from  tropical  forests 
may  only  become  apparent  after  many  years  of  observation. 

There  were  large  interannual  differences  in  apparent  forest-level  quantum 
efficiency  (a)  estimated  from  the  NEEday  data.  This  may  indicate  large  adjustments  in 
forest  structure  and  physiology  in  response  to  the  climatic  variation  among  our  study 
years.  Waring  et  al.  (1995)  concluded  that  seasonal  differences  in  both  LAI  and  canopy- 
level  quantum  efficiency  largely  controlled  productivity  from  a  deciduous  forest  in 
northeastern  US.  In  our  study,  the  estimated  a  for  2000  was  significantly  higher  than 
those  from  the  preceding  two  years  and  approached  the  theoretical  maximum  for  C3 
leaves  (Lloyd  et  al.  1995,  Farquhar  et  al.  1980).  However  eLAI  did  not  follow  this 
pattern  and  increased  only  1999  and  into  the  dry  season  of  2000.  That  neither  changes  in 
LAI  nor  a  were  related  to  annual  NEE  suggests  that  there  are  interactions  between  NEE, 
canopy  leaf  dynamics  and  climate  for  this  complex  tropical  wet  forest. 

The  relationship  between  NEEday  and  VPD  could  be  a  result  of  either 
physiological  or  physical  effects.  A  physiological  effect  could  be  stomatal  closure  in 
response  to  a  hydrologic  limitation,  either  high  VPDs  or  indirectly,  decreases  in  soil 
moisture  availability.  The  location  of  the  hydraulic  limitation  in  the  La  Selva  forest  is  not 
known.  A  physical  effect  could  be  through  though  modification  of  canopy  architecture 
through  premature  leaf  drop,  leaf  folding  or  changes  in  leaf  orientation.  Whole  forest 
canopies  do  not  frilly  saturate  even  at  frill  insolation  (Ruimy  et  al.  1995,  Wang  and 
Polglase  1995).  Changes  in  leaf  angle  or  leaf  closure  in  the  upper  canopy  allows 
penetration  of  light  to  deeper  canopy  layers,  allowing  for  increased  carbon  gain  in  the 
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lower  canopy.  This  offsets  the  effect  of  leaf  closure  or  changes  in  orientation  at  leaf  scale 
in  terms  of  light  response  at  the  ecosystem  level.  Only  8%  amount  of  the  time,  however, 
were  VPD  values  >  1  kPa  and  during  97%  of  the  daytime  (when  net  radiation  was  >  40  w 
m'^),  latent  heat  fluxes  were  greater  than  sensible  heat  fluxes  (P  <1,  unpublished  data). 
This  strongly  suggests  that  the  La  Selva  canopy  had  access  to  abundant  soil  water.  The 
only  exception  was  in  the  1998  dry  season,  when  30%  of  daytime  VPDs  were  >  1  kPa, 
precipitation  was  the  lowest  ever  recorded  (68,  38,  126  mm  monthly  total  rainfall  in 
January,  February  and  March,  respectively),  and  daily  mean  Ta  was  ~1  °C  above  the  long- 
term  average. 

The  1998  dry  season  was  at  the  end  of  the  1997-1998  warm-phase  El  Nino 
Southern  Oscillation  (ENSO)  and  was  warmer  and  drier.  During  December  1 998,  a  cold- 
phase  (la  Nina)  ENSO  brought  greater  precipitation,  cooler  temperatures,  and  lower  mean 
daily  insolation  (and  PPFD),  with  several  days  receiving  <  5  MJ  d"'.  Overall,  1998  was 
warmer  and  drier  during  the  dry  season,  but  had  more  precipitation,  cooler  temperatures 
and  reduced  light  during  the  latter  part  of  the  year  (compared  to  the  other  two  years,  2-4). 
The  greater  eLAI  in  1998  coupled  with  lower  a,  a  greater  portion  of  time  with  VPDs  >  1 
kPa  in  the  dry  season,  and  overall  lower  mean  daily  insolation  (13.3  MJ  d"',  2-4),  likely 
contributed  to  the  La  Selva  forest  being  close  to  carbon  neutral  in  1998. 
During  1999,  the  daily  insolation  was  well  above  the  long  term  trend,  but  in  November 
and  December  the  insolation  was  well  below  the  long-term  average  due  to  a  prolonged 
temporal,  suggesting  that  the  effects  of  increased  annual  incident  radiation  outweighed 
those  of  reduced  eLAI  and  a  prolonged  temporal,  making  this  forest  a  moderate  sink  of 
carbon. 
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Despite  these  seasonal  variations  in  climate,  I  did  not  find  any  seasonal  effects  on 
NEE  for  any  year.  Even  though  seasonal  displacements  of  the  ITCZ  alter  Hadley  cell 
circulations,  changes  in  individual  environmental  factors  do  not  necessarily  occur  in 
concert  as  a  result.  The  initial  passing  of  the  ITCZ  can  be  intermittent  and  there  can  be 
multiple  'false  starts'  (Hastenrath  1991).  Moreover,  the  northern  most  progression  of  the 
ITCZ  is  just  north  of  Costa  Rica  and  with  erratic  movement,  prolonged  periods  of  dry 
weather  can  occur  during  the  otherwise  wet  season  (Sanford  et  al.  1994).  Other  regional 
anomalies  can  occur,  as  during  1998  and  1999,  when  heavy  rains  persisted  fi^om 
December  into  January,  even  after  the  ITCZ  passed.  In  1998,  this  was  brought  about  by  a 
cold-phase  ENSO  event,  and  in  1999  by  a  prolonged  temporal. 

The  climatic  trends  observed  during  the  97/98  ENSO  were  tjpical  for  this  region 
(Waylen  et  al.  1 996b,  Cavazos  and  Hastenrath  1 990),  as  were  the  conditions  observed 
during  temporales  (Sanford  et  al.  1994).  The  97/98  ENSO,  however,  brought  the  highest 
temperatures  in  the  1 9-year  La  Selva  record.  The  apparent  increasing  frequency  and 
magnitude  of  ENSO  events  (Timmermann  et  al.  1999,  Corti  et  al.  1999),  may  have 
implications  for  interpreting  the  effects  of  climate  change  on  NEE,  as  at  La  Selva,  1 998 
had  the  lowest  estimated  NEE.  This  supports  the  findings  of  Tian  et  al.  (1998)  who 
concluded  in  a  modelling  study  that  variations  in  NEE  of  tropical  forests  are  controlled, 
in  part,  by  macro-level  changes  in  climate,  which  in  turn  are  driven  by  the  timing, 
frequency  and  magnitude  of  ENSO  events. 

A  possible  alternative  explanation  for  the  large  interannual  variation  in  NEE  at  La 
Selva  is  that  a  substantial  fraction  of  forest  is  in  an  early  successional  stage.  There  is  a 
high  frequency  of  treefall  gaps  in  this  forest,  even  though  La  Selva  is  below  most 


25 

hurricane  pathlines.  If  I  assume  that  the  forest  is  aggrading  carbon  during  its  stand  half- 
life,  which  was  estimated  as  77.3  y  (Lieberman  et  al.  1990),  and  has  gap  formation  rate  of 
0.96%  area  y"'  (Denslow  and  Hartshorn  1994),  then  74%  of  the  land  area  is  under 
recovery  at  any  time.  Moreover,  the  mortality  rate  for  1969-1982  was  2.03%  (Lieberman 
et  al  1990),  but  increased  to  4.77%  in  1997-1998  (Clark  et  al.  unpublished  data)  in 
upland  plots  associated  with  gaps,  suggesting  that  gap  formation  rate  also  increased  under 
climatic  conditions  imposed  by  both  ENSO  phases. 

Comparisons  with  Other  Tropical  Sites 

The  pattern  and  magnitude  of  NEE  at  La  Selva  were  similar  to  those  estimated  by 
eddy  covariance  at  three  tropical  moist  forest  sites  in  the  Brazilian  Amazon  (Reserva 
Ducke,  north  central  Amazon,  Fan  et  al.  1990;  Jani,  south  central  Amazon,  Grace  et  al. 
1995b;  Cuieiras,  north  central  Amazon,  Mahli  et  al.  1998).  Mean  maximum  daytime 
NEE  estimated  for  these  three  forests  ca,  -18-to-20  fxmol  CO2  m'^  s"'  and  mean  nighttime 
NEE  -5-7  ^imol  CO2  m'^  s''  (Table  2-5).  This  suggests  that  these  ecosystems  may  have 
similar  controlling  factors  on  NEE,  even  though  there  are  marked  differences  in  stand 
characteristics.  Mahli  et  al.  (1998)  hypothesized  that  cloudiness  on  insolation  were  the 
strong  determinate  of  NEEday,  which  may  explain  the  low  NEE  in  1998  at  La  Selva.  One 
difference  is  that  VPD  was  found  to  play  a  stronger  role  in  regulating  carbon  gain  at  the 
other  sites  than  at  La  Selva  (e.g.,  Mahli  et  al.  1 998).  These  other  three  sites  receive  < 
2500  mm  y"'  of  annual  precipitation  and  have  lower  soil  water  availability  (Hodnett  et  al. 
1996,  Tomasella  and  Hodnett  1996),  so  that  there  are  likely  to  be  greater  hydrologic 
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constraints  on  daytime  NEE  (Mahli  et  al.  1 998,  Williams  et  al.  1 998).  The  range  in 
annual  quantum  yield  at  La  Selva  overlaps  with  the  values  for  the  other  three  sites. 

Reduction  in  leaf  carbon  gain  can  occur  when  leaves  are  wet,  lowering  rates  of 
CO2  diffusion  into  the  leaf  by  a  factor  of  10"*  (Jones  1992),  thus  reducing  photosynthesis 
(Ishabashi  and  Terashima  1995).  Smith  and  McClean  (1989)  also  found  that 
photosynthesis  was  significantly  reduced  in  wet  leaves  that  had  a  wettable  cuticle,  but 
increased  on  leaves  that  had  non-wettable  cuticles.  The  ratio  of  species  with 
wettable  :non- wettable  leaves  is  at  La  Selva  not  known.  However,  because  eddy 
covariance  data  collected  during  rain  events  at  La  Selva  were  eliminated,  NEEday  may 
have  been  over-estimated  as  a  result.  La  Selva  experiences  rain  for  18%  of  the  time 
annually,  as  compared  to  only  13%  and  8%  at  Jani  and  Ducke,  respectively  (data  fi-om 
1994-96,  httpW.www.abracos.com).  For  this  reason,  the  likelihood  of  a  sampling  bias  are 
potentially  greater  at  La  Selva  due  to  a  larger  fi-action  of  data  removed  in  screening. 

The  range  of  our  estimates  of  annual  gross  ecosystem  production  (GEP,  Table  2- 
6)  overlaps  the  GPP  estimate  at  Cuieiras.  High  ratios  of  NEE:GEP  are  thought  to  exist 
only  with  forests  rapidly  accumulating  carbon  (e.g.,  0.02  for  an  old-growth  Pseudotsuga 
menziesii  forest,  compared  to  0.29  for  rapidly  a  growing  pine  stand.  Waring  and 
Schlesinger  1 985).  The  NEE:GEP  ratios  were  0. 1 8  and  0. 1 94  for  La  Selva  and  Cuieiras, 
respectively,  during  2000  and  1 995  indicating  that  these  systems  can  be  very  productive 
compared  to  other  forests.  All  of  the  tropical  sites  cited  here  are  considered  "old- 
growth",  and  assumed  in  the  past  to  be  near  steady  state  with  respect  to  carbon. 
However,  eddy  flux  data  fi-om  none  of  them  support  this  view  and  suggest  instead  that 
tropical  forests  may  be  net  sinks  from  0  to  much  as  -6 1  C  ha'*  y"'  (averaging  ca.  -2  t  C 
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ha''  y'').  However,  the  systematic  errors  discussed  earlier  collectively  would  tend  to 
decrease  the  amount  of  accumulated  NEE  into  the  forest,  since  the  carbon  would  be 
measured  when  it  enters  the  system  at  the  top  but  not  when  it  leaves  the  forest  (e.g., 
downslope  drainage  of  CO2,  loss  of  C  as  unmeasured  volatDe  organic  compounds,  etc.). 
Where  within  the  La  Selva  ecosystem  this  carbon  is  apparently  accumulating  remains  an 
open  question. 
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Table  2-1 .  Estimates  of  aerodynamic  parameters,  zero-plane  displacement  (d)  from  eq. 
2-1 ,  aerodynamic  roughness  length  (zo)  from  eq.  2-2,  and  u*  estimates  according  to 
stability  class  (L),  eq.  2-3,  from  La  Selva  Biological  Station,  zo,  d  and  u*  estimates  are 
median  values  ±95%  CI,  L  values  are  means  ±1  SE,  and  n  is  number  of  30-min  periods. 


Stability  class 

zo(ni) 

d(m) 

u*  (m  s"') 

L(m) 

n 

Unstable 

2.41  ±1.04 

21.5  ±1.80 

0.37  ±0.46 

-716  ±428 

57 

Slightly  unstable 

3.62  ±0.95 

19.4  ±1.91 

0.34  ±0.13 

-22  ±1.5 

51 

Neutral 

0.45  ±0.01 

22.1  ±6.73 

0.11  ±0.05 

1  ±0.6 

84 

Stable 

0.44  ±0.21 

23.0  ±1.00 

0.22  ±0.05 

814  ±553 

186 
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Table  2-3.  Annual  and  seasonal  differences  in  estimated  leaf  area  index  (eLAI) 
m'^  from  La  Selva,  Costa  Rica  (S.F.  Oberbauer  unpublished  data). 


Season  Year 

Mean  (median)  ±  1  SE 

Dry  1998 

3.85  (3.96)  ±0.19 

Wet  1998 

4.85  (4.79)  ±0.11 

Dry  1999 

2.71  (2.52)  ±0.13 

Wet  1999 

3.76  (3.84)  ±0.11 

Dry  2000 

3.48  (3.51)  ±0.07 

Wet  2000 

3.43  (3.30)  ±0.13 
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Table  2-6.  Between  year  measures  of  productivity  and  ecosystem  efficiency  from  La 
Selva,  Costa  Rica  and  Cuieiras,  Brazil.  Data  from  La  Selva  are  from  eq.  2-5  across  VPD 
classes  and  a  fixed  nighttime  NEE,  as  noted  in  the  text.  Units  of  GEP  and  NEE  are 
expressed  as  MT  C  ha''  y''. 


site 

La  Selva,  Costa  Rica 

Cuieiras,  Brazil 

year 

1998 

1999 

2000 

1995 

GEP 

28.41 

30.6 

33.9 

30.4 

NEE 

-0.09 

-1.66 

-6.1 

-5.9 

NEE:GEP 

0.004 

0.055 

0.18 

0.194 

Note:  The  Cuieriras  dataset  is  from  Mahli  et  al.  1998,  and  Malhi  and  Grace  2000.  Data 
from  La  Selva  is  based  on  NEEmodeied. 
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Figure  2-1 .  The  relationship  between  the  normalized  power  spectra,  for  the  three  wind 
vectors  against  the  wavenumber.  Data  are  averaged  from  73  30-min  periods  from 
January  10  to  Februray  28th,  1999.  A)  Mid-day,  unstable  conditions,  beginning  at  1200- 
1230  p.m,  and  B)  Mid-night,  stable  conditions,  beginning  at  2300-2330  p.m.  Error  bars 
are  +/-  1  SE. 
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Figure  2-2.  Diurnal  time  series  of  the  spectra-based  correction  factor, 

for  CO2  fluxes  from  La  Selva  as  calculated  from  Eq  2-5.  Data  are  90-min  running  means 

centered  on  the  30-min  interval.  Data  are  from  1998-2000.  Error  bars  are  +/-  1  SE. 
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Figure  2-3.  Relationships  between  u*  and  daytime  NEE,  (A),  and  between  B)  nighttime 
eddy  covariance  or  storage  fluxes,  and  u*  (B).  Values  of  nighttime  NEE  were  averaged 
across  intervals  of  0.025  m  s"',  except  for  the  righthand  most  eddy  covariance  point, 
which  was  averaged  from  all  data  between  0.4  <  u*  <  0.54  m  s'\  Error  bars  are  +/- 1 SE. 
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Figure  2-4.  Diurnal  characteristics  of  A)  storage  fluxes,  and  B)  line-averaged 
temperature  and  non-rotated  vertical  windspeed  from  La  Selva.  All  data  are  from  1 998- 
2000.  A  90-niin  running  average  was  used  with  each  estimate  centered  on  the  30-min 
interval  as  indicated.  Sample  size  for  storage  and  line-averaged  temperature  were  5566 
and  for  vertical  windspeed  16333.  Error  bars  are  +/-  1  SE. 
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Figure  2-5.  The  mean  diumal  pattern  of  above-canopy  eddy  co variance  (EC)  and  storage 
(below  canopy)  flux,  NEE  (the  sum  of  the  fluxes)  and  PPFD  from  La  Selva  (1998-2000). 
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Figure  2-6.  Main  effects  of  environmental  variables  on  daytime  NEE, 

where.  A)  is  the  light  response  of  NEE  for  ail  data  (1998-2000),  and  B)  and  C)  show  the 

residuals  from  the  light  response  function  in  relation  to  temperature  and  VPD, 

respectively. 
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Figure  2-7.  NEE  as  a  function  of  PPFD  across  a  gradient  of  vapor  pressure  deficits 
for  years  1998,  1999  and  2000  fi-om  La  Selva,  Costa  Rica. 
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Figure  2-7.  continued. 
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Figure  2-7.  continued. 
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Figure  2-8.  The  relationship  between  NEEnight  and  line-averaged  (below-canopy) 
temperature.  NEEnight  estimates  are  the  aggregated  averages  from  each  nighttime  30-min 
period  from  the  whole  study  period  (Figure  2-4A).  Temperature  is  an  average  from  six 
measurements  through  the  canopy  profile  from  27  to  0.5  m. 
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Figure  2-9.  Cumulative  NEE  from  La  Selva,  Costa  Rica  for  1998-2000. 
Estimates  were  calculated  using  A)  the  results  from  the  light  response  equation  across 
VPD  classes  and  a  fixed  estimate  for  nighttime  respiration,  and  B)  direct  estimates  of 
daytime  and  nighttime  NEE  with  gaps  filled  using  data  from  NEEmodeled 


CHAPTER  3 

ENERGY  BALANCE  AND  MODELED  EVAPOTRANSPIRATION  FOR  A  WET 
TROPICAL  FOREST  IN  COSTA  RICA 

The  objective  of  this  chapter  was  to  examine  the  annual  and  seasonal  effects  of 

atmospheric  environment  (albedo,  net  radiation,  vapor  pressure  deficit)  and  the  surface 

controls  on  the  energy  balance  of  a  wet  tropical  forest  in  Costa  Rica. 

Methods 

Meteorological  Data 

All  measurements  referred  to  here  were  collected  fi-om  September  1997  to 
December  2000.  Instrumentation  for  measuring  air  temperature,  relative  humidity,  bulk 
precipitation,  and  net  radiation  were  mounted  at  the  top  of  a  42  m  tower  (Upright,  Inc., 
Selma,  CA).  Prior  to  March  1  1999,  air  temperature  (Ta)  was  measured  with  a  CS500 
probe  (Campbell  Scientific,  Inc.,  Logan  UT)  installed  within  a  radiation  shield,  and 
linearly  back-corrected  to  fit  the  response  of  the  aspirated  temperature  sensor  (R  =  0.98). 
Ta  was  measured  with  platinum  resistance  temperature  detector  (100  Q  platinum  RTD, 
Omega  Engineering,  Stamford,  CT)  mounted  in  an  aspirated  shield.  Relative  humidity 
was  also  measured  with  the  CS500  probe,  and  rainfall  with  a  tipping  bucket  rain  gauge 
(model  TE525,  metric,  Texas  Electronics,  Dallas  TX).  Net  radiation,  R„,  was  measured 
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from  March  1,  1999  to  December  2000,  with  a  closed-ceU  thermopile-style  sensor  (NR- 
Ute,  Kipp  and  Zonen,  Delft,  the  Netherlands).  From  August  1,  1997  to  March  1, 1999  R„ 
was  measured  with  a  Fritschen-style  sensor  (model  Q-7.1,  Radiation  Energy  Balance 
Systems,  Seattle,  WA).  All  data  collected  with  the  Q-7.1  were  linearly  back-corrected  to 
fit  the  response  of  the  NR-Ute  (R^  =  0.97)  and  corrected  for  advected  sensible  heat. 
Another  net  radiometer  (model  CNR.l,  Kipp  and  Zonen,  Delft,  the  Netherlands)  was  also 
used  to  estimate  albedo  during  February-April,  July,  and  September  2000.  Two 
resistance  grid  type  leaf  wetness  sensors  (model  237-L,  CampbeU  Scientific,  Inc.)  were 
mounted  at  26  and  2  m,  and  histograms  of  relative  wetness  were  compiled.  To  mimic  the 
wetting  of  leaves,  the  sensors  were  coated  with  a  layer  of  flat,  off-white  latex  paint. 

Soil  heat  flux  plates  (model  HFT-3,  Radiation  Energy  Balance  Systems)  were 
installed  at  a  depth  of  5  cm,  in  each  of  three  1  m  x  1  m  plots  >  20  m  distance  apart  near 
the  base  of  the  tower.  Atmospheric  pressure  was  measured  at  ~3  m  (PB105,  Vaisala, 
Helsinki,  Finland).  All  meteorological  data  were  collected  at  5  sec  intervals  and 
compiled  as  30-min  averages  with  a  datalogger  (CRIOX,  Campbell  Scientific  Inc., 
Logan,  UT).  Instruments  were  cleaned,  leveled  as  necessary,  and  recalibrated  according 
to  manufacturers'  instructions. 


Energy  Balance  Estimates 

An  ecosystem-level  energy  balance  can  be  estimated  by 
R„=AE  +  H  +  G  Eq.  3-1 

assuming  horizontal  homogeneity,  where  Rn  is  net  radiation,  XE  is  the  latent  energy  flux, 
H  is  the  sensible  heat  flux,  and  G  is  the  soil  heat  flux  (all  units  are  W  m'^).  Both  H  and 
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XE  were  estimated  by  the  summation  of  both  eddy  co variance  (above-canopy)  and  below- 
canopy  fluxes,  expressed  by, 


where  w',  q',  and  9'  are  the  deviations  of  instantaneous  values  from  a  running  mean  of 
vertical  windspeed  (m  s"'),  specific  humidity  (mmol  mol''),  and  virtual  temperature  (°C), 
respectively,  and  Zx  is  measurement  height  (m).  w'q'  and  w'0'  are  the  turbulent 
exchanges  of  water  vapor  and  heat  as  estimated  by  eddy  covariance  method,  and  the 
second  term  in  each  equation  is  the  storage  flux  in  the  air  column.  The  storage  fluxes  are 
noted  as  X,Estor  and  Hstor,  respectively.  The  convention  used  here  is  that  negative  values 
correspond  to  a  flux  into  the  forest  from  the  atmosphere. 

The  eddy  covariance  system  was  comprised  of  a  sonic  anemometer  (K-probe, 
Applied  Technologies,  Inc.,  Boulder,  CO),an  infra-red  gas  analyzer  (IRGA,  model  Li- 
6262,  Li-Cor,  Linclon,  NE),  -60  m  of  tubing  (4.8  mm  ID  Teflon)  with  the  inlet  co- 
located  with  the  sonic  probe,  a  laptop  computer,  and  a  pump  to  pull  air  through  the  tubing 
at  ~  8  1pm.  The  sonic  anemometer  measured  the  wind  velocities  in  three  dimensions  at 
10  Hz,  where  w  is  vertical  windspeed,  and  u  and  v  are  the  two  horizontal  windspeed 
components.  The  anemometer  was  also  used  to  estimate  temperature,  9s,  as  a  function  of 
the  speed  of  sound  and  changes  in  air  density  (excluding  water  vapor,  i.e.,  "virtual 
temperature").  The  IRGA  measured  the  concentration  of  water  vapor  at  ~8  Hz.  The 
laptop  and  a  flux  software  program  from  McMillen  (1988),  with  a  fixed  lagtime  (14.3  s), 
were  used  to  collect  raw  eddy  covariance  data  files.  Covariances,  wind  and  scalar 


Eq.  3-2 


Eq.  3-3 
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statistics,  and  coordinate  rotations  were  calculated  in  real  time  at  10  Hz.  Because 
transport  of  mass  and  energy  by  turbulence  occurs  between  frequencies  of  1-10  Hz, 
within  the  inertial  sub-range  (Kaimal  and  Finnigan  1994),  and  because  the  IRGA's 
response  time  is  ~8  Hz,  fast-Fourier  analyses  were  applied  to  correct  for  frequency  loss. 
Protocols  for  accuracy,  precision,  quality  control  and  assurance  were  used  as  defined  by 
the  AmeriFlux  Science  Plan  (http://cdiac.esd.oml.gov/programs/ameriflux/scifhtm). 

The  storage  fluxes  were  estimated  by  adding  together  the  changes  of  heat  through 
the  forest  profile  in  both  the  air  column  below  42  m  and  in  the  foliage  (Eq.  3-3).  Water 
vapor  was  sampled  from  6  inlets  at  0.5,  7.3,  1 1.95,  16.55,  21.2,  and  27.6  m  on  the  tower. 
Solenoids  switched  the  flow  (~3  1pm)  from  each  inlet  through  a  second  IRGA  for  5  min 
during  each  30-min  period.  For  H,  temperature  through  the  profile.  Op,  was  measured 
with  platinum  RTDs  housed  in  radiation  shields  and  co-located  with  each  inlet;  when 
sampling  occurred,  the  airflow  acted  to  aspirate  the  platinum  RTDs,  approximating  true 
8. 

Because  leaves  have  small  thermal  inertia  but  significant  amounts  of  H2O, 
changes  in  leaf  H  were  estimated  by: 


where  Hieaf  is  the  below-canopy  leaf  heat  flux  (W  m"^),  S^vt  is  the  specific  weight  of  leaf 


one-sided  leaf  area  based  on  LAI  from  polyculture  plantations  at  La  Selva  (4.5  m  m" ,  S. 
Bigelow  and  J.Ewel,  pers.  comm.).  It  was  assumed  that  the  profiles  of  Hstor  and  XEstor 
were  similar  throughout  the  flux  source  area.  G  was  estimated  as  a  30-min  average  of  the 


0.5 


Eq.  3-4 


water  (kg  H2O  m'^  leaf  area),  Cpw  is  the  heat  capacity  of  water  (J  kg'  K' ),  and  Larea  is  the 
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3  soil  heat  flux  plates.  The  contribution  to  G  by  the  change  in  soil  heat  storage  was 
ignored  because  soil  temperature  varied  little  in  a  30-min  period,  soil  thermal  properties 
for  these  soils  were  not  known,  and  the  change  in  soil  heat  storage  was  expected  to  be 
negligible. 

Evapotranspiration  Model 

Evapotranspiration  estimates  were  partitioned  into  whole  forest  transpiration  and 

the  evaporation  of  intercepted  precipitation.  XEj  was  calculated  using  the  Penman- 

Monteith  equation; 

M„+PaC[e,iTa)-ea(T^)]ga 
XE  =   Eq.  3-5 

/l[A  +  Kl  +  ")] 
8b 

where  XE  is  latent  energy  flux  (W  m"^),  A  is  the  is  the  rate  of  increase  in  saturated  water 
vapor  pressure  with  temperature  (kPa  K""),  pa  is  the  density  of  air  (kg  m"^),  Cs  is  the 
saturated  water  vapor  pressure  at  Ta,  ea  is  the  ambient  water  vapor  pressure  (kPa),  ga  is 
the  aerodynamic  conductance  (mol  m'^  s"'),  A,  is  the  latent  heat  of  vaporization  (J  kg"'),  y 
is  the  psychometric  constant  at  25  °C  (0.0665  kPa  K''),  gb  is  the  bulk  canopy  conductance 
(mol  m"^  s'').  To  change  units  of  energy  to  depth,  XE  was  multiplied  by  a  conversion 
factor  that  included  molar  volume  (mol  m"^)  and  weight  (kg  mol'').  Evapotranspiration 
depth  is  noted  in  the  text  as  Ej. 

A  positive  momentum  flux  into  the  canopy  was  assumed  hence,  aerodynamic 
conductance  was  estimated  by, 
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k^u  Eq.  3-6 
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where  k  is  Von  Karmen's  constant  (0.40),  d  is  the  zeroplane  displacement  (m),  Zm  is  the 
aerodynamic  roughness  length,  and      and  Th  are  the  diabatic  correction  factors  (m)  for 
momentum  and  sensible  heat,  respectively  (Yasuda  1988,  Arya  1988).  Zeroplane 
displacement  and  aerodynamic  roughness  changed  with  stability  and  were  empirically 
estimated  for  this  study  period  (Table  2-1,  Loescher  et  al.  submitted).  Diabatic  correction 
factors  are  a  function  of  stability,  where  in  stable  conditions, 

=^>„  =61n(l+0  ^'l-^-^ 
and  in  unstable  conditions, 

^  -21n[— =  0.6T^ 

with  C|  is  a  stability  parameter  ratio  of  convective  to  mechanical  turbulent  production  , 

Z42j^  Eq.  3-9 

L 

where  L  is  the  Monin-Obukov  length  (Equation  2-3). 
Bulk  canopy  conductance,  gb  was  estimated  by, 
_     gg        Ag?-  Eq.  3-10 

CpAE 

where  D42  is  the  specific  humidity  deficit  at  measurement  height  (kg  kg"'). 

When  gaps  in  the  measured  XE  occurred,  both  ga  and  gb  were  empirically 
modelled  by  using  relationships  with  horizontal  windspeed  (for  ga),  and  the  vapor 
pressure  deficit  (VPD)  and  Rn,  (for  gb,  Martin  et  al.  1997,  Wright  et  al.  1996).  gb  was 
then  normalized  to  unity,  and  the  upper  limit  to  VPD  determined  (Jarvis  1976,  Livingston 
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and  Black  1987).  This  limit  function  was  used  to  estimate  a  theoretical  maximum,  gmax, 
by  increasing  gb  as  though  VPD  was  not  limiting.  Then  gmax,  in  turn,  was  related  to  R„, 
assuming  that  maximum  conductance  would  take  place  with  0  VPD  and  high  R„.  A 
dimensionless  decoupling  coefficient,  Q,  was  used  to  determine  the  relative  effects  of  ga 
and  gb  on  evapotranspiration  (Jarvis  and  McNaughton  1986); 

A/ +2 

Q  =   Eq.3-11 

/r  Sb 

Evaporation  of  free  water  in  the  canopy  was  modeled  using  a  Rutter-type  model 
(Calder  et  al.  1986).  Because  canopy  water  storage  increases  exponentially  to  a 
maximum,  with  increases  in  precipitation  from  a  single  rain  event,  Ei  was  modeled  every 
30-min  by, 


t  =  0 


Eq.  3-12 

where  t  is  min,  C  is  canopy  water  depth.  When  the  canopy  was  wet,  Ex  was  estimated 
with  gb  set  to  0.  C  was  estimated  by, 

C^C,^A\-e    '^'''P)  Eq.3-13 


where  Ccap  is  the  maximum  canopy  capacity,  mm,  k  is  a  unitless  canopy  fill  coefficient , 
SP  is  the  cumulative  amount  of  water  that  fell  during  a  30-min  interval.  An  empirical 
estimate  of  1 .53  mm  was  used  for  Ccap  (Loescher  et  al.  2002).  The  stemflow  component 
of  interception  was  ignored  because  it  was  assumed  to  be  a  small  volumetric  flux,  i.e.,  < 
2%  of  rain,  (Schroth  et  al.  1 999,  Neal  et  al.  1993).  A  value  of  0.28  for  k  from  a  broadleaf 
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plantation  forest  at  La  Selva  was  used  (Bigelow  2001).  If  modeUed  Et  was  greater  than 
the  remainder  of  canopy  free  water,  the  remaining  depth  was  included  as  transpiration, 
Et.  Evapotranspiration  was  then  the  sum  of  Ei  and  Et. 

A  second  method  of  estimating  XEr,  the  Priestly-Taylor  equation,  was  used  to 
compare  with  the  Penman-Monteith  results.  The  Priestly-Taylor  equation  (Eq.  3-14) 
simplifies  the  transfer  process  that  is  explicit  in  Eq.  3-5,  and  in  doing  so,  is  thought  to  be 
appropriate  for  large-scale,  well-watered  vegetative  canopies,  like  those  typically  found 
in  the  wet  tropics  (Priestly  and  Taylor  1972). 

Ai:^=aR„[-^]  Eq.3-14 

where  a  is  a  coefiBcient  estimated  by  fitting  the  model  results  to  empirical  measures  of 
A.E  from  Eq.  3-2.  Monteith  (1981)  estimated  an  average  a  on  a  theoretical  basis  as  1 .26, 
but  values  observed  over  rough  canopies  have  varied  greatly  (Jones  1992). 

Results 

Above-Canopy  Environment 

The  Monin-Obukov  (M-0)  stability  length  did  not  differ  with  year  or  season.  On 
a  diurnal  basis,  M-0  length  was  neutral  (~  0  m)  during  the  night,  and  decreased  during 
the  day  time  until  1400  when  the  boundary  layer  became  weakly  unstable  (M-0  ~  -125  m 
Figure  3-1).  After  1400  h,  the  M-0  length  sharply  increased,  and  the  boundary-layer 
became  weakly  stable  (-100  m)  at  1600  h,  but  then  returned  to  neutral  conditions  by 
nightfall.  Friction  velocity  (u*)  was  -0.1  m  s"'  during  the  night  and  increased  during  the 
morning  hours  with  convective  turbulence,  u*  decreased  after  solar  noon  and  continued 
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to  decrease  throughout  the  afternoon  due  to  the  dissipation  of  turbulent  kinetic  energy 
(Figure  3-1). 

Total  daily  Rn  ranged  from  1 .47  to  27.54  MJ  d"'  during  the  measurement  period, 
and  differed  among  years  (using  a  general  linear  model  with  an  a  =  0.05,  p  <  0.0001, 
Figure  3-2)  with  mean  daily  totals  for  1998-2000  of  13.31  ±0.028,  17.48  ±0.050,  and 
15.33  ±0.040  MJ  d''  (mean  ±1  SE),  respectively.  Rn  also  varied  significantly  with  season 
(p  <  0.0001).  Mid-day  albedo  did  not  change  seasonally,  and  ranged  from  0.118-0.135  of 
incident  short  wave  radiation. 

Enei^  Balance 

Diurnal  temperature  and  water  vapor  profiles  (Figure  3-3)  followed  trends  similar 
to  those  of  other  forests  (Shaw  et  al.  1988).  Heating  of  the  air  column  during  the  day  was 
greater  with  height,  i.e.,  there  was  a  positive  temperature  gradient.  However,  counter 
gradients  were  often  observable  between  21  and  27  m,  where  the  leaf  area  was 
concentrated.  Cooling  during  the  night  often  produced  neutral  or  slightly  negative 
gradients,  often  with  warmer  temperatures  at  ground  levels.  Negative  or  neutral  water 
vapor  gradients  were  observed  all  times,  with  counter  gradients  present  during  non-rain 
days  between  1 1  and  21  m  height. 

Soil  heat  fluxes  followed  very  similar  diurnal  patterns  throughout  the  year,  and 
ranged  between  ±  1 6  W  m~^  at  any  point  in  time,  with  negative  flux  into  the  system 
during  the  daytime  (Figure  3-4A).  At  night,  XEstor  was  ~  3-5  W  m'^.  A  larger  A,Estor 
efflux  occurred  in  the  early  morning  hours,  presumably  from  convective  winds  mixing 
the  below-canopy  airspace  and  evaporating  free  water.  X,Estor  flux  continued  to  be 
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positive  throughout  the  afternoon,  but  was  more  variable,  with  mean  daytime  values 
ranging  from  ~6  to  7  W  m"^  This  flux  remained  positive  throughout  the  day  across  all 
years.  Mean  nighttime  storage  of  sensible  heat  (Hstor  +  Hieaf)  ranged  from  ~  1-3  W  m" , 
and  became  negative  in  the  early  hours  as  the  air  space  increased  in  temperature.  The 
maximum  average  Hstor  +  Hieaf  was  ca.  -18  W  m"^,  which  occurred  at  -0800  h  when  the 
convective  boundary  layer  was  developing.  Hstor  +  Hieaf  increased  and  became  positive  at 
-1400  h,  which  coincided  with  a  late  afternoon  weakly  stable/unstable  boundary  layer 
(Figure  3-1).  The  below-canopy  environment  continued  to  lose  sensible  heat  until  -1900 
h,  when  neutral  canopy  conditions  prevailed.  The  total  storage  flux  (Figure  3-4B)  was 
similar  to  that  of  Hstor  +  Hieaf,  but  daytime  fluxes  were  ameliorated  by  release  of  water  to 
the  atmosphere.  During  nighttime  neutral  conditions,  the  flux  was  -8-10  Wm".  The 
daytime  minima  was  -  -5  W  m'^  and  the  peak  efflux  was  17  W  m"^,  which  occurred 
during  weakly  stable  conditions  (Figure  3-4B). 

The  average  30-min  XE  was  greater  than  H  for  all  daytime  hours  and  across 
seasons  and  years  (i.e.,  H/X.E  =  p  <  1 .0,  Figure  3-5A).  A  linear  model  that  included 
second-order  effects  of  year,  season,  VPD  and  R,,,  explained  79%  of  the  total  variation  in 
H  +XE  (Table  3-1).  Because  Rn  and  VPD  are  auto-correlated  and  VPD  did  not  explain 
any  additional  variation,  it  was  removed  from  the  linear  model.  Rn  alone  accounted  for 
69%  and  68%  of  the  variation  in  H  and  XE,  respectively. 

I  could  only  close  the  energy  balance  (Eq.  3-1)  to  within  32-50%,  with  the 
exception  of  periods  with  a  wet  canopy  during  the  dry  season  of  1998.  In  general,  H  + 
XE  +  G  underestimated  R,,,  particularly  when  Rn  values  were  <  400  W  m'^  (Figure  3-6). 
The  most  likely  reason  for  this  was  that  XE  estimates,  which  also  contributed  to  greater 
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variance  in  p  estimate  before  1000  h  (Figure  3-5 A).  For  this  reason,  XE  was  empirically 
estimated  by  >.E  =  Rn  -  H  for  model  comparisons,  rather  than  the  A,E  estimate  in  Eq  3-2 
(e.g.,  Twine  et  al.  2000).  Bowen  ratio  calculated  by  H/(Rn-H)  were  -50%  smaller  than 
those  estimated  by  IE  in  Eq  3-2,  and  the  largest  contribution  of  H  occurring  between  0900 
and  1200  h  (Figure  3-5B),  coinciding  closely  with  the  observed  diurnal  pattern  of  Tg. 

Modelled  Conductance  and  Evaporation 

ga,  was  linearly  correlated  with  horizontal  windspeed  (Figure  3-7A).  Values  of 
\|/H  and  vj/M  were  small  and  all  other  parameters  were  relatively  constant.  We  expected 
the  upper  boundary  of  normalized  gb  (Figure  3-7B)  to  be  negatively  related  to  VPD, 
because  of  the  negative  physiological  response  of  canopies  to  VPD  (Martin  et  al.  1997), 
which  was  found  for  VPDs  >  0.5  kPa.  gmax  was  hyper  bo  lically  related  to  Rn  (Figure  3- 
7C).  Estimated  values  of  gb  were  always  less  than  ga  (Figure  3-8A).  At  dawn,  gb  was  ~1 
mol  m"^  s'*  imtil  0800  h,  then  steadily  decreased  during  the  day,  approaching  0. 1  mol  m"^ 
s"'  at  dusk.  The  minimum  ga  was  1 .5  mol  m"^  s'',  which  increased  to  >  2  mol  m'^  s'' 
during  mid-day  (1000-1500  h).  Modeled  values  of  ga  and  gb  behaved  similarly  to  those 
derived  using  Eq.  3-2  (Figure  3-8B).  Using  conductances  in  the  Penman-Monteith 
equation  during  periods  when  the  canopy  was  dry,  ga  explained  44%  of  the  variation  in 
^E  (Figure  3-9).  Q  ranged  from  -0.3  during  the  night  to  -0.7  by  0830  h.  For  the 
majority  of  daytime  hours  (0600-1600),  Q  was  >  0.5  (Figure  3-10). 

The  Penman-Monteith  equation  explained  68%  of  the  observed  variation  in  Rn  - 
H,  but  overestimated  Rn  -  H  by  -28%  (Figure  3-1 1  A).  In  contrast,  the  more  simple 
Priestly-Taylor  relationship  accounted  for  98%  of  the  observed  variation  (Figure  3-1  IB). 
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Annual  ETft  ranged  from  1892  mm  in  1998  to  2294  mm  in  1999,  and  from  54%  to  66% 
of  bulk  precipitation  (Table  3-2).  Mean  daily  ETpt  rates  were  also  lower  in  1998  and 
greatest  in  1999  (5.19  to  6.29  mm  d'',  respectively).  Interception  loss  was  greatest  in 
2000,  with  an  annual  total  of  708  mm,  accounting  for  18%  of  bulk  precipitation. 

Discussion 

Ecosystem  Enei^  Dynamics 

Daytime  profiles  of  temperature  and  water  vapor  in  the  upper  canopy  often 
showed  counter  gradients  that  were  the  result  of  winds  that  did  not  frilly  penetrate  the 
entire  canopy.  The  higher  canopy  vegetation,  between  21-27  m,  was  a  physical  barrier  to 
transpiration  from  below.  As  a  result,  the  apparent  A,Estor  at  night  was  due  adiabatic 
cooling  and  consequent  condensation.  Storage  fluxes  can  contribute  substantially  to  the 
overall  ecosystem  energy  balance  when  R„  is  small  or  during  the  night,  but  on  a  diurnal 
basis  this  was  a  very  a  small  component  (~2  %). 

Bowen  ratios  (P)  consistently  <  1  indicated  that  water  was  not  limiting  XE  at  any 
time  during  the  year.  Soils  in  the  wet  tropics  generally  do  not  exert  hydraulic  limitations 
on  XE  (De  Bruin  1983).  This  is  likely  also  the  case  at  La  Selva,  where  soils  have  high 
water-holding  capacity  and  high  hydraulic  conductivity  (Weitz  et  al.  1 997,  SoUins  et  al. 
1994).  Pentaclethra  macroloba,  the  dominant  tree  (42%  of  the  basal  area),  closes  its 
stomata  and  leaves  in  the  late  afternoon  (~  1530  h).  I  could  not  detect  associated  changes 
in  either  H  or  XE  due  to  changes  in  P.  macroloba  physiology,  suggesting  that  either 
single  species  control  on  energy  partitioning  or  cannot  be  detected  at  the  ecosystem-level 
for  diverse  tropical  wet  forests,  the  eddy-covariance  technique  is  insensitive  to  species- 
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level  responses  in  these  forests ,  or  that  the  time  of  day  that  P.  macroloba  closes  its 
leaves  is  inconsequential  to  energy  balance. 

The  lack  of  energy  balance  closure  at  this  site  is  probably  attributable  to  both 
losses  of  high  and  low  frequency  signals.  This  loss  is  more  apparent  when  R,,  <  400  W 
m"l  A  combination  of  the  long  sampling  tube  combined  with  high  humidity  conditions 
may  have  led  to  some  in-tube  mixing.  Also,  the  scale  lengths  for  H  and  IE  were  highly 
variable  and  may  not  have  been  fully  captured  using  eddy  covariance  (Mahli  1 996). 

Similar  degrees  of  non-closure  have  been  found  at  other  tower  flux  sites  with 
similarly  large  aerodynamic  roughness  (~  2  m),  such  as  an  old-grow4h  pine  stand  in 
Oregon  (to  within  20-30  %,  Anthoni  et  al.  2000),  and  an  old-growth  conifer  stand  in 
Washington  ( to  within  10-35  %,  M.  Falk  and  J.  Chen,  pers.  comm.).  Studies  from  other 
tropical  have  not  reported  the  degree  of  energy  closure  for  periods  greater  than  1  day,  but 
it  is  likely  that  similar  results  would  also  be  obtained  for  these  also  complex  natural 
forests.  A  high  degree  of  closiire  for  30-min  intervals  seems  to  require  aerodynamic 
roughness  lengths  less  than  -0.5  m  to  minimize  viscous  effects,  less  than  saturated  air  for 
long  periods,  and  smaller  sensor  separation  distances  (e.g.,  Gholz  and  Clark,  in  press). 

Conductances  and  Other  Limits  to  Annual  Energy  Fluxes 

The  general  diurnal  patterns  found  for  ga  and  gb  are  similar  to  those  for  other 
tropical  forests  (Bigelow  2001,  Wright  et  al.  1996,  Shuttleworth  1984).  Using  the  same 
calculations,  our  gb  estimates  were  almost  identical  to  those  found  by  Wright  et  al.  (1996) 
for  a  tropical  forest  in  Brazil,  suggesting  similar  physical  and  physiological  controls  on 
from  these  two  neotropical  forest  canopies.  Bigelow  (2001)  examined  gb  for  three 
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monocultural  plantations  at  La  Selva  and  found  that  all  the  species  had  higher  late 
afternoon  rates  (-0.5-0.8  mol  m"^  s"')  than  found  for  the  old-growth  forest.  However  two 

2     1  • 

species,  Cedrela  odorata  and  Cordia  alliodora,  had  late  afternoon  gb  <  0.2  mol  m'  s"  in 
the  dry  season.  The  lower  afternoon  gb  that  we  observed  was  likely  an  integrated 
response  of  P.  macroloba  closing  its  stomata  and  folding  its  leaves  closed  by  1630  h  or 
later.  Significant  differences  were  observed  in  gb  between  years  and  seasons,  which 
followed  the  same  trends  as  Rn.  Between  year  differences  in  Rn  explained  much  of  the 
interannual  variability  observed  in  measured  and  modelled  XE.  Increased  rainfall  in  2000 
increased  the  absolute  amount  of  interception,  but  not  the  fi-action  of  rainfall  intercepted. 

The  direct  effects  of  gb  on     could  not  be  determined  because  these  variables 
could  not  be  independently  measured.  However,  low  values  of  Q  (e.g.,  <  0.3)  during 
night  and  early  morning  and  evening  hours  indicate  that  these  periods  are  the  only  times 
when  physiological  control  over  XE  occurs,  likely  due  to  the  opening  and  closing  of 
stomata  and  leaves.  Values  of  Q  >  0.5  suggest  that  mid-day  X,E  is  controlled  more  by  Rn 
and  ga  than  gb.  ga  explained  44  %  of  the  variation  in  XE  during  times  when  both 
conductances  were  used  to  calculate  ETpm  (i.e.,  dry  canopy  conditions).  The  importance 
of  ga  in  controlling  ETpm  increased  fiirther  because  of  the  very  high  precipitation  at  La 
Selva.  For  example,  32  %  of  the  time,  the  upper  canopy  was  wet. 

The  Priestly-Taylor  relationship  for  well- watered  conditions  described  the  ET 
dynamics  of  this  tropical  forest  quite  well.  The  fi-action  of  available  energy  used  for  ETpj 
was  similar  fi-om  year-to-year  suggesting  a  thermodynamic  constraint  on  H  that  limits 
maximum  daily  virtual  temperature  (Calder  1986).  Wright  et  al.  (1996)  and  Calder  et  al. 
(1986)  also  found  that  accounted  for  a  large  fi-action  of  Rn  (>0.80)  from  a  humid 
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Amazonian  and  Javanese  forest,  respectively,  as  did  Bigelow  (2001)  for  the  monocultural 
plantations  at  La  Selva  (0.79-0.90). 

In  a  previous  study  in  the  old-grovvth  forest  at  La  Selva,  Luvall  (1984)  determined 
that  the  energy  required  for  evapotranspiration  exceeded  Rn  by  25%.  Although  this 
seems  counterintuitive,  it  may  very  well  be  true.  This  phenomenon  has  been  observed 
over  crops  (Ham  and  Hiehnan  1991)  as  well  as  other  tropical  forests  (Jones  1992, 
Shuttleworth  1989,  Calder  1986).  The  likely  explanation  is  that  energy  is  locally 
advected  into  the  flux  field.  This  is  certainly  possible  at  La  Selva,  where  mean  daytime 
wind  direction  is  -90  and  the  fetch  is  ~2  km,  beyond  which  the  landscape  is  dominated 
by  pastures,  crops  and  patches  of  secondary  forests  extending  for  ~  60  km  to  the 
Caribbean  shore.  Often,  small  (<  1  km  in  width)  afternoon  convective  cells  deposit  rain 
heterogeneously  across  the  landscape.  Advection  of  drier  air  masses  with  greater 
evaporative  demand  is  possible,  particularly  during  Luvall's  study  in  the  early  80's  when 
much  of  the  land  in  the  Costa  Rican  coastal  plain  was  being  converted  from  forests  to 
agriculture.  Estimates  of  ETpr  reported  here  did  not  exceed  available  Rn,  suggesting 
advection  was  not  significant,  in  this  case,  although  we  cannot  rule  out  the  possibility  that 
advection  contributed  additional  energy  towards  the  overall  balance. 

The  magnitude  of  canopy  capacitance  is  in  large  part  a  fiinction  of  physical 
surface  area  of  a  canopy  (Waring  and  Schleshinger  1985).  At  La  Selva,  high  epiphytic 
loads,  bromeliad  tanks  and  arboreal  soil  mats  can  contribute  capacitance  and  may  not 
have  been  fiilly  accoimted  for  in  our  estimates.  We  used  a  fixed  estimate  of  capacitance 
of  1.53  mm  The  fraction  of  rainfall  that  was  intercepted  seemed  consistent  from  1998  to 
2000,  and  also  with  Luvall's  study,  which  suggests  that  the  canopy  surface  area  at  La 
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Selva  is  saturated  for  a  large  portion  of  time,  that  the  relative  annual  amount  of 
throughfaU  is  constant,  and  that  surface  area  does  not  change  to  any  appreciable  degree 
over  time. 
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Table  3-1.  First  order  regression  parameters  for  the  energy  balance  closure  from  1998- 
2000  across  seasons  and  canopy  conditions  from  a  wet,  tropical,  old-growth  forest  at  La 
Selva,  Costa  Rica,  y  =  Yo  +  b*x,  y  is  X,E  +  H,  x  is  Rn,  Yo  and  b  values  are  mean  ±1SD,  n 
is  number  of  30-min  averages,  and  all  regressions  were  significant  at  the  p  <  0.0001  level 

unless  otherwise  noted.  ^_ 

year  Yo  b  n  R 


1998 


Dry  canopy/dry  season 

12.0  ±3.8  p 

=  0.0016 

0.52  ±0.02 

429 

0.72 

Wet  canopy/dry  season 

38.0  ±6.2 

1.07  ±0.05 

362 

0.59 

Dry  canopy/wet  season 

12.7  ±2.6 

0.66  ±0.01 

2182 

0.73 

Wet  canopy/  Wet  season 

7.9  ±2.3  p= 

0.0006 

0.63  ±0.11 

1867 

0.66 

1999 

Dry  canopy/dry  season 

16.0  ±2.3 

0.63  ±0.01 

2372 

0.83 

Wet  canopy/dry  season 

19.4  ±2.5 

0.57  ±0.01 

1362 

0.75 

Dry  canopy/wet  season 

13.1  ±2.3 

0.55  ±0.01 

2356 

0.74 

Wet  canopy/  Wet  season 

22.6  ±5.0 

0.50  ±0.02 

483 

0.73 

2000 

Dry  canopy/dry  season 

7.3  ±4  p  =  0.07 

0.61  ±0.01 

1220 

0.81 

Wet  canopy/dry  season 

17.5  ±5.5  p 

=  0.0014 

0.6  ±0.02 

256 

0.84 

Dry  canopy/wet  season 

17.5  ±2.6 

0.57  ±0.01 

1868 

0.68 

Wet  canopy/  Wet  season 

23.1  ±4.1 

0.68  ±0.2 

872 

0.71 
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Table  3-3.  Annual  fraction  of  time  that  the  canopy  was  wet  at  two  heights  in  the  La  Selva 
canopy. 

year  Height  of  leaf  wetness  sensor 


25  m  2  m 

1998  042  059 

1999  0.17  0.57 

2000  0.23  0.48 
All  years            0.32  0.57 
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Time  (hour) 


Figure  3-1 .  Diumal  relationships  of  Mction  velocity  and  Monin-Obukov  length 
over  an  old-growth  tropical  forest.  Data  are  averages  using  all  data  from  1998-2000. 
Intervals  A,  B,  C,  and  D  are  neutral,  unstable,  weakly  unstable,  weakly  stable  boundary 
layers,  respectively.  Error  bars  are  +/-  1  SE. 
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Figure  3-2.  Cumulative  net  radiation  for  1998-2000  over  an  old-growth  forest 

in  La  Selva,  Costa  Rica.  Daily  means  were  derived  using  first-order  regression.  Each 

year  was  significantly  different  at  the  a  =  0.05  level,  p  <  0.001,  and     >  0.99. 
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Figure  3-3.  Typical  diurnal  changes  in  below-canopy  temperature  and  water  vapor 
profiles  fi-om  an  old-growth  forest,  La  Selva,  Costa  Rica.  Data  are  median  values  for  aU 
of  1999.  Error  bars  are  +/- 1  SE. 
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Figure  3-4.  Diumal  patterns  of  storage  energy  fluxes  from  an  old-growth  wet  tropical 
forest  in  Costa  Rica.  Data  are  mean  values  from  1998-2000  with  error  bars  +/-  1  SE. 


70 


UJ 

£^ 
X 

o 

(0 

i_ 

c 
a> 

o 
m 


1.8 
1.6 
1.4 

1.2  ^ 

1.0 
0.8 
0.6 
0.4 

0.2 

0.0 
-0.2 
-0.4 


o 


c 

9) 
O 

m 


If 


2000 

1      1  1999 
1      1  1998 

..ill 

Mil 

Mill 

B 

fittl. 

1                 J                 1                 1  1 

600 


800 


1000 


1200 


Time  (hour) 


1400 


1600 


Figure  3-5.  Daytime  Bowen  ratios  for  each  year.  Estimates  of     were  calculated  by  A) 
equation  3-2,  and  B)  Rn  -  H,  with  H  determined  by  equation  3-3.  Data  are  measured 
median  (+  1  SE)  values  with  net  radiation  >  40  W  m'^. 
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Figure  3-6.  The  relationship  between  net  radiation  and  estimated  energy  flux 
for  1998-2000.  Estimated  energy  flux  includes  modelled  contribution  from  the 
understory. 
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Figure  3-6.  continued 
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Figure  3-6.  continued 
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g,  =  1.37  +  0.66'u 
=  0.91 


Figure  3-7.  Emprical  relations  to  model  both  aerodynamic  and  bulk  conductance, 
independant  variables  were  averaged  at  different  intervals,  i.e.,  0.25  m  s' , 
0.025  kPa,  and  50  w  m-2,  for  A,  B,  and  C,  respectively.  All  data  were  median 
values  wdth  +/- 1  SE. 
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Figure  3-8.  The  diumal  relationship  of  aerodynamic  and  bulk  conductance  calculated  by 
A)  collected  eddy  covariance  data,  and  B)  modelled  based  on  Figure  7.  Data  are  median 
values  with  +/- 1  SE. 


Figure  3-9.  The  relationship  between  aerodynamic  conductance  and  latent  energy  flux 
from  La  Selva  during  1999.  Data  are  empirical  estimates  from  Eq.  3-6  (ga),  and  from  Eq. 
3-2  (XE). 
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Figure  3-10.  Diurnal  changes  in  the  decoupling  coeflBcient,  Q,  from  an  old-growth  wet 
tropical  forest,  Costa  Rica.  Data  is  from  1997-2000.  SE  are  typically  <  0.006. 
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Figure  3-11.  Relationship  between  empirical  and  modeled  estimates  of  the  latent  energy, 
where  A)  is  modeled  using  a  Penman-Monteith  equation,  and  B)  using  Priestly-Talyor 
equation  with  a  dry  canopy,  and  C)  using  a  Priestly- Taylor  equation  with  a  wet  canopy. 
All  graphs  used  data  from  1998-2000.  An  a  of  1 .24  was  found  for  graphs  B  and  C.  All 
slopes  were  significant,  p  <  0.0001 .  y- intercepts  had  p- values  of  0.5,  0.04  and  0.04  for  A, 
B  and  C,  respectively.. 


CHAPTER  4 
CONCLUSIONS 

Diurnal  patterns  of  NEE  at  La  Selva  followed  trends  similar  to  those  observed 
elsewhere,  with  turbulent  fluxes  dominating  in  the  daytime,  and  large  storage  fluxes 
contributing  largely  in  the  early  morning.  Daytime  NEE  was  a  function  of  both  abiotic 
(PPFD,  VPD  and  temperature)  and  biotic  (quantum  efficiency  and  eLAI)  factors.  VPD 
limited  NEEday  when  values  were  above  >  1  kPa,  but  this  only  occurred  over  a  small 
percentage  of  time.  There  was  a  positive  correlation  between  nighttime  respiration  and 
temperature  based  on  diurnal  averages  of  all  NEEnight  measurements.  We  used  a  fixed 
NEEnight  value,  to  estimate  annual  NEE.  There  was  a  big  difiference  between  our  two 
alternative  estimates  of  annual  NEE.  But  regardless  of  the  calculation  method  used,  the 
results  indicate  that  there  is  a  large  interannual  variation  in  NEE  at  La  Selva,  related  to 
large  scale  regional  climate  dynamics. 

Energy  balance  closure  of  this  lowland  tropical  rainforest  was  ~  74%.  Storage 
fluxes  contributed  very  little  (-2%)  to  the  overall  daily  energy  balance.  Conductances 
foUowed  similar  trends  to  those  found  of  other  tropical  forests.  Daytime  XE  was  almost 
always  greater  than  sensible  heat,  suggesting  that  the  trees  in  this  forest  are  in  contact 
with  sufBcient  ground-water  reserves  to  minimize  hydraulic  stress.  Rn  was  a  large 
determinant  for  the  annual  energy  flux,  suggesting  that  the  Priestly-Taylor  model  for 
evapotranspiration  is  more  appropriate  in  the  tropics  than  the  Penman-Monteith  model. 
The  general  rule  holds  true  that  Rn  ~  XE. 
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